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THE CONCEPTS OF NOTIONAL SHIP AND NOTIONAL VALUE 
IN LOGISTIC CAPABILITY STUDIES INVOLVING MERCHANT SHIPS* + 


Ralph B. Hunt and Erling F. Rosholdt 
The George Washington University 
Logistics Research Project 


INTRODUCTION 

A primary logistic planning operation is that of determining how many ships, or freight 
cars, or trucks, or planes, or even buses would be required to transport a given quantity of 
cargo (dry, liquid, or human) from one point to another according to a particular lift require- 
ment schedule. It is generally impossible or at the least impracticable to relate the capacity, 
speed, and cargo-handling characteristics of each specific type of carrier to the demands of a 
lift requirement schedule. For studies involving ships, a way around the problem has been 
found through the use of the concepts of a "notional" ship and of the "notional" value of an actu- 
al ship. This paper describes in some detail these two concepts with the hope that a similar 


approach may be found useful with other kinds of cargo carriers. 


THE CONCEPT OF A NOTIONAL SHIP 

The cargo capacity, speed, and cargo-handling characteristics of merchant ships are 
quite diverse. Within certain classes of ships of similar design, e.g., the Liberty ships, 
there is a similarity in the characteristics; but, considering all existing merchant ships, it 
would be nearly a hopeless task to include each ship's actual characteristics in the parameters 
of most logistical planning problems. Accordingly, the concept of a "notional" ship seeks to 
provide a way to handle the variety of ships' characteristics in a feasible manner. This con- 
cept (which is not original with the authors but the origin of which is not known) conceives of a 
reference ship of a standard cargo capacity, speed, and cargo-handling capability in terms of 
Which the actual ships would be expressed. This comparison is expressed as the "notional 
value" of an actual ship and is a dimensionless ratio. In its usual sense, the "notional value" 
implies that the characteristics of an actual ship give it more or less "value" than the refer- 
ence (notional) ship. Certain factors in this concept are discussed below. 

Consider two ships which operate over the same route. With respect only to the 


Capacities of the ships, other things being equal, that ship having the greater cargo capacity 


*Manuscript received June 25, 1959. 
tResearch supported by the Office of Naval Research under Task NR 047-001. 
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will be able to carry more payload per trip and therefore will have greater "value" to its 
owners. Considering now only speed, that ship with the higher speed will complete the trip in 
the shorter time and will be available more quickly for another payload. In other words, its 
"turnover" rate is higher and, therefore, the resource (the ship) has greater "value." Hence, 
its "value" is a function of its speed. The ability to complete a trip quickly and have a high 
"turnover" value depends not only on speed but also on the ability to spend a minimum time in 


port, since a ship earns nothing while alongside a dock. Assuming port facilities are adequate, 


the time in port that is spent handling cargo will depend, at least to a first approximation, on 
the capacity of the ship (quantity to be handled) and the number of hatches available. (A few 


ships may have extra large hatches so that more than one gang of stevedores can work a hatch; 


also, more than the usual number of cargo winches may be available. However, the productiv- 
ity of stevedores and longshoremen are possible offsetting quantities to the hatch size, and only 
very few ships have additional winches. Hence, we submit that, in the general case, it is rea- 
sonable to assume that the loading time is inversely proportional to the number of hatches. ) 
Therefore, the ship with the more rapid cargo-handling capability will spend less time in port 
and therefore complete its trip more quickly, making possible a higher "turnover" rate and 
increasing its ''value.'' Perhaps "value" should be considered synonymous with "usefulness" 
in our frame of reference. 


The effects of speed and port time appear in the ship's turnaround time (Eq. (1)). 


(1) TA = Pp + Dg + Pt + Dg 


where (all time in days) 
TA = turnaround time, 
= time spent in port of origin, 
= one-way steaming time, 
= time spent in terminus or destination port. 
Notional value is then expressed in the form 
c, (TA) 
(2) NV, = Cc. ° (FA) 
n 
where subscripts a = an actual ship, 
n = a notional ship, 
r = denotes an actual route (2-way), 


(?) denotes a notional ship over an actual route (2-way), 


(2) denotes an actual ship over an actual route (2-way), 


(TA) = ship turnaround time in days, 
NV = notional value, 


C = capacity. 
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Equation (2) simply says that the notional value of an actual ship will vary directly as 
the cargo capacities and inversely as the turnaround times, i.e., greater capacity and lower 
turnaround time in an actual ship will increase its ''value" over the notional ship. 


For actual computation purposes, Eq. (2) must be expanded to the form given below: 





Dy 
a 24Sn 
NVa = ro . 
n D, 


+ 
24S, 


(dry cargo or bulk liquid), 
where (all time in days) 


subscripts n notional ship 
a actual ship 


r actual route 


notional value (freighter or tanker), 

cargo unloading time, 

cargo loading time, 

(W, + Wy + Up + Ly), where Wo waiting time at port of origin, 


Wt waiting time at port of 
destination, 


and if ballast is required, 


Up = ballast unloading time, 
Lp = ballast loading time, 
Dry two-way distance in nautical miles, 
24 hours/day, 
s speed in knots, 
D,/24S = two-way steaming time in days, 
c capacity (deadweight tons (DWT) for dry cargo and liquid cargo) 
H number of hatches (applies only to dry cargo ships). 
Note: Omit the correction terms [ (C,/C,) * (H),/H,) | for liquid-cargo ships since, for 
planning purposes, the capacity has only a minor effect on the pumping time when (U) and (L) 
are expressed in days. Increase both (U) and (L) by one day for "monster" tankers. 
Two things in Eq. (3) require further discussion. First, the equation gives the no- 


tional value of an actual ship for a specific value of the distance, D,- The point to be em- 





phasized is that the notional value of an actual ship is a function of the distance (route) over 





Which it operates. If port times and capacities remain the same, the notional value will vary 


with the distance. Second, the unit of capacity is given as the deadweight ton (DWT = 2240 lIb.). 


The unit here is a matter of convenience. Any other convenient set of units for cargo capacity 
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would be equally satisfactory so long as both capacities are expressed in the same units. In 


naval usage, the measurement ton! is commonly used. For a specific set of conditions the 


M/T equivalent of the DWT capacity may be varied according to the judgment of the planner. 
For example, if C, is fixed as 10,800 DWT, the equivalent in M/T may be taken as 7,500, 
10,000, etc., depending on whether the conditions of the plan will permit normal cargo hand- 
ling or require special handling (e. g., priority unloading). Then the requirements, which 
are given in M/T, can be converted to notional ship equivalents, which would then be based on 
M/T units. An equivalence between barrels and tanker capacity in DWT is used in a similar 
manner. 


For human cargo (passengers), Eq. (3) takes the form of 


Dr 
Ca (TA)n Ca | U, + Li + Ky + ae 


NV s = — 
. C n (TA da Ch Ca Dy 





r (U, + La) so K, + 
| Cy 24S, 


where the capacity units are numbers of persons and the other symbols have the same meaning 
as in Eq. (3). 

The procedures for determining the number of notional ships required to meet a given 
lift-requirement schedule have been described in detail* and will not be discussed here. 
After the number of notional ships required over each route has been determined, it is then 


necessary to convert this figure to the number of actual ships required. 


"NOTIONAL" VALUE AND "NOTIONAL" DISTANCE 

The total actual ship requirements are obtained by applying the concept of the "notional" 
value of an actual ship which has been discussed earlier. It must be remembered that the 
notional ship requirements available at this stage are for specific routes. As has been 
pointed out, the notional value of an actual ship is a function of the distance (route) over 
which it operates. There are two ways of handling this factor. One is to determine the 
notional value of each actual ship for every route for which the number of notional ships 
needed were determined. The other is to decide on a standard reference distance ("notional" 
distance) and convert the notional ship requirements, as calculated for the different routes, 
to the notional ship equivalent at the notional distance selected. Since by using a different 


distance we are changing the turnaround time of the notional ship, then the ratio of the 


1 A measurement ton is the space available for cargo measured in units of 40 cu ft to the long 
ton (2240 lb). A capacity of 10,000 M/T is the same as 400, 000 cu ft measured to the inside 
of the cargo battens, on the frames, and to the under side of the beams. 

R. B. Hunt, and E. F. Rosholdt, ''Time-Phasing Relationships in Determining Merchant 
Ship Requirements for Overseas Cargo Lift: Part 1. The Export-No Pipeline System, '' The 
George Washington University Logistics Research Project, Washington, D. C., Serial 
T-101/59, 15 May 1959. 
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number of the notional ships required for the notional distance to the number required for the 


individual route distance will be 
(TA) p+ 


(TA 
{ ng 





(5) ( oes) r 248 


(Qy P+ Dnd 
24S 
[24S - P+ Dy} (Q)r 


(Qnd = “948. P+ Dna 





where subscripts notional ship 
= notional distance (2-way) 
actual route (2-way) 


notional ship over actual route 


notional ship over notional distance 


number of notional ships 

turnaround time, in days 

port time, in days (same for both turnarounds) 

2-way distance (nautical miles) 

notional speed, knots 

24 = hours/days 

After converting the number of notional ships required over each route to the equivalent num- 
ber required over the notional distance, then the notional value of the actual ships available 
can be determined according to Eq. (2). 

The use of a "notional'’ distance is most convenient when it is necessary to make a 
capability study of "notional" ships required over all routes. By converting the notional ships 
required over each route to the equivalent number required over the notional distance, the 
sum of the notional ships based on the notional distance can be obtained and compared with the 
cumulative notional value of available actual ships. Table 1 illustrates this comparison for a 
hypothetical example. 

TABLE 1 


Effect of Applying Notional-Distance Concept to Hypothetical Notional-Ship 
Requirements over Different Routes Requiring Different Turnaround Times 





Actual Notional Turnaround Turnaround Time | Equivalent Notional 
Two-Way | Ships Required at | Time for Actual} Over Notional Ships at Notional 
Distance Actual Distance Distance Distance Distance of 

(n.m.) (Q),. (days) (8000 n. m. )(days) 8000 n.m.*(Q)ng 





4000 100 iy ee : 60.6 





6000 100 22.7 ; 80.3 





9000 100 31.0 ; 109 .9 




















_ 11000 100 36.6 . 129 .6 














~Q, 400 380.4 











*Determined by Eq. (5). 
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Thus, a total notional value (at 8000 miles) of 380.4 ships (not 400) must be satisfied 


by available actual ships with a cumulative notional value (at 8000 miles) which is equal to 


or greater than this sum. If the notional value(at 8000 miles) of a particular type of actual 


ship was, say, 1.3, then the number of actual ships of that type that would be required 
would be 380. 4/1.3 = 293. 
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A POSTSCRIPT TO "DYNAMIC PROBLEMS IN THE THEORY OF THE FIRM"'* 


Harvey M. Wagner 


Stanford University 


INTRODUCTION 

In a previous paper | 3] the combined use of traditional economic constructs and 
dynamic programming was suggested for solutions to several intertemporal problems in the 
theory of the firm (that produces a single commodity. )} Asa specific application of these 
techniques, a method was presented for solving a dynamic version of the economic lot size 
model |4]- The purpose of this postscript is: (i) to demonstrate that the algorithm suggested 
for the dynamic economic lot size model is also applicable to situations in which the cost 
curves (which may differ from period to period) have nonincreasing marginal costs as a 
function of output. 2 e.g., because of quantity discounts, and (ii) to describe a forward 
algorithm for solving optimal pricing and output problems in which the cost curves have 
nondecreasing marginal costs as a function of output .° As in [3], the approach emphasizing 
the use of analytical techniques familiar to economists is continued. ‘ Throughout we assume 
that each marginal revenue relation is a nonincreasing function of the corresponding period 
sales and that initial and ending inventory are zero. 

We consider finding a sequence of output qo(t) and prices, which in turn determine 
sales gg (t), that maximizes total profits over the entirety of periods t=1, 2, .-., T. The 
period marginal revenue and cost functions are denoted by MR [qs(t)] and MC [q,(t)], and 
the unit cost of carrying an item of inventory from period t to period t+1 is denoted by it. 
We postulate that at each period t, MR [qg(t)] < M (finite); that for qs(t) sufficiently large, 
MR [ q,(t)] < 0; and MC{[qo(t)] >0 for all qo(t). 

It will be helpful to recall the previous system of diagrams [3, p.57], arraying each 
period's marginal revenue and cost curves, constructed such that the second period's curves 
are shifted downward by the amount ij, the third period's by the amount ij+ ig, ..., and the 


ph period's by the amount ij + ig+ ... + ip. ,. All geometric arguments throughout the paper 


*Manuscript received January 13, 1959. 
1 The reader should recognize the parallel existing between ''production" and ''ordering for 
inventory;'' for ease of exposition, we frame our discussion in the language of "production". 
2Figuratively, "increasing returns to scale." 

3Figuratively, "decreasing returns to scale." 

There are differences of opinion as to whether or not such geometric devices are more 
awkward, theoretically or computationally, than a recursive equation formulation [1]. Our 
predilection in this paper in motivated by a desire to bridge the gap between traditional eco- 
nomic analysis and dynamic programming. 
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will subsume that the marginal curves are shifted appropriately. We permit the marginal 
functions to have horizontal segments (representing a constant unit cost or price over a 
range of values for output or sales) and adopt the convention of connecting two "adjacent" 
horizontal segments by a vertical segment.” It should be noted that in the commonly con- 
sidered case of a fixed amount demanded Ss for each period t (at some fixed price p,) the 


marginal revenue curve is represented as a step function 


Pt if Gg (t) < St 


MR [ag (t)] -{ 


0 otherwise. 


In case of a capacity constraint on period production, the marginal cost curve is depicted as 


rising vertically at this point. 


THE CASE OF NONINCREASING MARGINAL COSTS 

Both [ 3 and 4] contain a discussion of the model in which the marginal cost curve is 
horizontal® and identical in all periods, and nonnegative setup costs are incurred if production 
takes place. The suggested forward algorithm is briefly described as: 

At period t* (starting with t* = 1) compute the profits from producing in period 
t** (t** = 1,2,...,t*), and filling all qg (t), t=t**, t** + 1,...,t*, by production in 
period t**. The method for determining q, (t) is illustrated in some detail in [ 3, pp. 
56-61]. In brief, if production in period t** is to be allocated to sales qg (t), t = t**, 
t** + 1,...,t*, then the curves MR [q, (t)] for the relevant periods are aggregated (by 
adding the relations horizontally), the total amount to be produced and sold is determined 
by the intersection of MC | qg (t**)| and the aggregate marginal revenue curve, say at the 
value MR, and this total quantity is distributed for sale among the periods according to the 
output associated with the value MR on each period's MR Ld, (t)] function. Add to the profit 
figure so determined for periods t**, t** + 1,...,t*, the profits using an optimal policy in 
periods t = 1,2,...,t** - 1 considered by themselves.’ Select producing in the t** that 
offers the maximum total profits; this yields an optimal policy for the first t* periods. 
Continue until t* = T. 

That this algorithm leads to an optimal solution is based on the fundamental proposi- 
tion [| 4, p. 91]: There exists an optimal program such that at any period t the firm need 


not both produce and enter the period with previous periods’ production. 


In other words, if a marginal curve is in part a step function, the ''rise'' in the step is 
drawn as well as the "flat." 

Consequently capacity constraints are ignored. 

Since the algorithm starts at t* = 1, these profit figures have been previously computed. 
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We assert that the same algorithm leads to an optimal solution if the marginal cost 
































curves are nonincreasing and not necessarily identical in all periods; as before, nonnegative 
setup costs may be included. The reason underlying the assertion is that the fundamental 
proposition continues to hold. We sketch a proof: 

Suppose a trial solution indicates in period t* that both inventory produced during pre- 
vious periods enters and production takes place. Specifically, let | P] denote the set of periods 
supplying inventory to period t*. Consider the current value of marginal cost for each t e| P] 


and add to this number the holding costs charged for carrying an additional unit of production 


BS ee cee eae te 


from t to t*. Let MC' be the minimum value of these sums and t' be the associated period. If 
MC' is less than the MC for the last unit produced in period t, we may revise the trial solution 
so as to increase production in period t' and eliminate it in period t. If MC’ is not less than 
the MC of the last unit produced in period t, we may revise the trial solution so as to increase 
production in period t equal to the amount of incoming inventory and eliminate the correspond- 
g ing production in [ P]. Such alterations do not increase the total cost over that of the trial 
S ’ solution because of our hypothesis of decreasing marginal cost functions. 
tion Be We should recognize that our nresent assumption about marginal costs is too weak to 
prove the Planning Horizon Theorem | 4], which in part states that if it is optimal to produce 
‘ : in t* when periods 1 through t* are considered by themselves, then it is sufficient to consider 
programs for all T periods such that production takes place in period t*. In Table 1 we pre- 
sent an example in which the proposition is violated; the optimal two-period program is to 
4 produce in period 2, whereas the optimal three-period program is to produce only in period 1. 
fi Consequently, in the general case of nonincreasing marginal costs, the addition of a new 


period t* + 1 may cause a major revision in an optimal plan for periods 1 through t*. 





























= TABLE 1 
fit 
i Period 1 Period2 | Period 3 
Setup cost 12 | 2 | 5 
Unit Carrying 0 0 | 0 
charge 
sll Unit cost 8 | 9 | 10 
Selling price 15 | 15 | 15 
| | | 
Amount demanded 0 3 | 20. | 
ee 
: THE CASE OF NONDECREASING MARGINAL costs.® 
=] In [3] we described a recursive method for solving the class of problems in which 
# marginal cost is non-decreasing. Essentially, the technique involved finding a first trial 
i. a 


This assumption is meant to preclude set-up costs. 
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solution for all T periods; if the plan proved to be infeasible,’ finding a second trial solution 
to a smaller number of periods; and so forth, until finding a feasible solution for periods 
t=1,2,..., t* < T. Ift*< T, the method is repeated for periods t= t* + 1,..., T. 

Here we note that the ingenious method of Johnson [2] extends immediately to our 
problem to provide an alternative mode of attack. Specifically, it is possible to construct an 
optimal sales and production pattern, utilizing an algorithm starting at t= 1, and succes- 
sively adding the data for each following period. 

We define a "cascade function, " to be superimposed on the system of diagrams arraying 
each period's marginal curves, as a step function having a single horizontal step for each 
period and such that the heights of the steps are nonincreasing over time. Figure 1 shows an 


example of a cascade function for a three-period model. 


$ we (1) $ $ 






































° s° 
Figure 1 - Cascade function (Cj, C2, C3, C4, C5) 

The algorithm for constructing an optimal solution proceeds as follows: 

Step 1. In period 1 define the height of the step as the value of the marginal cost curve 
at the output for which period-1 marginal cost and revenue curves intersect. If the intersec- 
tion is at a vertical segment of the MC curve, then the height is defined as the value of MC 
on the immediate left of the intersection point.!0 

The provisional values for qs (1) and go (1) are determined as the output at the inter- 
section of the MR and MC curves. If the intersection is nonunique, then we assume that 
the largest salable output on the locus of intersection is selected. 

Step 2. The height of the cascade step for period 2 is provisionally set at the height of 
the step for period 1. The trial value for qo (2) is the smallest output such that MC [ao (2)| 
equals the height of the trial step, and for qs (2) it is the largest sales figure such that 
MR [as (2) | equals the height of the trial step. 


Case (i). If qo (2) = qs (2), no alteration is made in the provisional values. 


Case (ii). If qo (2)>qs(2), then the height of the step for period 2 must be lowered until 


equality obtains between the amounts sold and produced in period 2. In lowering the step, if 





9i.e., at some period, cumulative sales exceeded cumulative production. 
107, may happen that at qo (1) = qg (1) = 0 the MC curve lies above the MR curve, in 

which case it is never profitable to sell any item in period 1; the height of the cascade step is 
then defined as the value of the MC curve on the immediate right of q, (1) - 0. 
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a horizontal section of the MR curve is reached, the sales figure should be increased as 
much as possible toward meeting the equality condition. If a horizontal section of the MC 
curve is reached, output should be curtailed as much as possible toward meeting the equality 
condition. We adopt the rule that if horizontal sections in the MR and MC functions are 
met simultaneously, the increase in sales takes precedence over the decrease in output. 
Case (iii). If q9(2) < qg(2), then in this instance cumulative output over the two 
periods is not sufficient to meet cumulative demand. If possible, increase output in period 
2 at the same marginal cost with a view toward meeting the condition that cumulative output 
equals cumulative sales. Should the condition remain unsatisfied, if possible increase output 
in period 1 at the same marginal cost. Should the condition then remain unsatisfied, if pos- 
sible decrease sales in period 2 at a value of marginal revenue equal to the provisional height 
of the step. Should the condition continue to be unsatisfied, if possible decrease sales in 
period 1 at a value of marginal revenue equal to the provisional height of the step. Should the 
condition still remain unsatisfied, then raise the height of the step, maintaining the cascade 
property by raising the height of the step for period 1 and altering mutatis mutandis q)(t), 


g(t), t= 1,2, until equality obtains between cumulative sales and production. 





In raising the height of the step, if horizontal segments inthe MR and MC curves 
are encountered, an increase in output (starting at the latest period being considered) takes 
precedence over a decrease in sales. 

Step t*. The height of the cascade step for period t* is provisionally set at the height 
of the step for period t* - 1. As in Step 2, trial values for q)(t*) and q,(t*) are deter- 
mined. 

Case (i). If qo(t*) = qg(t*), then no alteration is made in the provisional values. 

Case (ii). If q )(t*) > ag(t*), then the alteration given in Step 2, Case (ii), is applied 
here. 

Case (iii). If qg(t*) < qg(t*), then if possible increase output in period t* at the 
Same marginal cost, with a view toward meeting the condition that cumulative output 
(t= 1,2,..., t*) equals cumulative sales. Should the condition remain unsatisfied, if possible 
increase output in periods t* - 1, t* - 2,..., at the same marginal cost. 11 Should the 
condition continue to be unsatisfied, if possible decrease sales in period t*, t*- 1, ...., 
ata value of marginal revenue equal to the provisional height of the step. Finally, should the 


condition still remain unsatisfied, then raise the height of the step, maintaining the cascade 





property by raising the height of the step for period t* - 1 and any previous periods as they 
become eligible and altering mutatis mutandis Go (t), g(t), t= 1,2, ..., t*, until equality 
obtains between cumulative sales and production. As in Step 2, Case (iii), increasing output 


takes precedence over decreasing sales whenever such simultaneous alternatives present 
themselves. 


| Sei . , ; ' 
1.€., for all immediately preceding periods having the cascade function at the same height 
as the current provisional value. 
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We outline an inductive proof for the validity of the algorithm. Suppose at period 
t* < T the ulgorithm has yielded an optimal feasible schedule for the first t* periods; this 
assumption holds for t* = 1. We now must verify that the procedure provides an optimal 
feasible schedule for periods 1 through t* + 1. If Case (i) or (ii) holds, the algorithm pro- 
vides a feasible schedule for periods 1 through t* + 1, for in essence period t* + 1 is 
considered separately; no revision in the previous periods' schedule is made, and in period 
t* + 1, qo(t* + 1) = qg(t* + 1) by construction. Since at this level of output and sales in 
period t* + 1 the marginal revenue and cost are no greater than the level of the cascade 
function in previous periods, there is no gain in profits to be made by carrying inventory into 
period t* + 1. 

If Case (iii) holds, q,(t* + 1) < qg(t* + 1) at the value of the cascade function for 
period t*. Therefore profits may be increased, even if the previous periods’ schedules are 
unaltered, by producing more and selling less than these trial amounts for period t* + 1. In 
making such revisions in output and sales, following Johnson | 2], we must also consider 
additions to production in previous periods and, similarly, reductions in previous periods 
sales. If the cascade function is shifted as indicated by the algorithm, we maintain the con- 
dition that there is no advantage of increasing further the level of inventory entering period 


t* + 1. Because of the assumptions concerning the shape and boundedness of the marginal 


revenue and cost curves, raising the cascade function eventually provides a feasible solution: 


at any period t < t* + 1, cumulative production increases and cumulative demand decreases; 
the upward shifitng eventually terminates, for a cascade value above M implies that 


cumulative demand is zero. 
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Performance assessment is a valuable tool for determining whether or 
not asystem is behaving satisfactorily and for isolating troubleareas if 
the performance is unsatisfactory. This type of evaluation is especially 
useful for systems that are in their initial stages of operation. If a new 
system represents a large-order combination of men and machines, 
development of a satisfactory method for evaluating such a system is 
not aneasy task. Experimental simulation involving actual (or laboratory) 
operation appears to furnish a feasible and realistic basis for such 
evaluations. By use of this approach, a model has been developed that 
seems to be acceptable for most evaluation purposes. This model is of 
a somewhat elementary mathematical form, being a linear sum ofthe con- 
tributions for important subdivisions with respect to system components 
and time; however, the model iselaborate in the sense that many inputs 
rely heavily on experience -judgment considerations. Explicitly, a 
measure of effectivenessis developed whichis based on a comparisonof 
what the system did with what it should have done, as determined by the 
recorded results of the simulation exercise. Anassessmentof a system's 
over-all behavior can be obtained by performing an evaluation for each 
ofa representative set of the important types of situations for this system. 











INTRODUCTION 

Many different evaluation viewpoints can be used for the assessment of a system. The 
viewpoint adopted here is that a system is behaving ideally if it accurately performs all of its 
intended functions at every given time of importance. Then the evaluation problem can be 
approached on the basis of the times that are considered to be important, a suitable breakdown 
of the system into basic components for each of these times, and the performance accuracies 
of the functions to be carried out by these components at these times. 

For new large-scale systems involving complicated combinations of men, machines, 
and perhaps electronics, interactions are virtually always very difficult to assess quantita- 
tively until operational experience is gained. Moreover, the number of interactions is usually 


so huge that their examination for all possibilities of interest is not feasible. Also, 
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many of the system activities may involve human actions, which are very difficult to predict 


(even on a probability distribution basis). Consequently, an experimental procedure in which 





































is used the actual personnel and equipment, combined with adequate recording of the inputs 
and outputs of the system during the experiment, appears to furnish as satisfactory a basis 
as any for a first-stage evaluation. A comprehensive recording of the experimental exercise 
can be used to determine what actually occurred and what should have occurred for the com- 
binations of times and basic components that are considered. 

Since the times chosen for consideration are seldom of anywhere near equal impor- 
tance and, for a given time, the same is true for the basic components, relative-worth 
functions are introduced for weighting the comparisons between actual and ideal behavior. 
These relative-worth functions are based on experience-judgment considerations and fulfill 
two purposes. One purpose is to assure that the viewpoint and model adopted are not | 
inconsistent with those for larger systems which include the system evaluated as a subsystem 
The other purpose is to make certain that the more important events receive the greater 1 
emphasis in the evaluation measure of effectiveness. 


The model adopted calls for computation of the worth of the system performance for 


each time considered and for the basic components considered at these times. The measure I 
of effectiveness is taken to be the accumulated worth of the system behavior for the exercise t 
being evaluated. That is, it is the sum of the worths for all the times of importance and cor- y 


responding basic components. An increase in the value of the measure of effectiveness indi- 
cates an improvement in system performance. j 
The stipulated additivity of worths for the subdivisions according to times and basic 
components has the advantage of simplifying the model but the disadvantage of tending to 
mask important interrelations among these subdivisions. Fortunately, this masking effect E 
can be at least partially eliminated by appropriate choice of the worth functions. The 
recorded results of an exercise should furnish an indication of the implications of inaccurate 
component behavior. A combination of these implications with experience and judgment 
considerations should yield worth-function expressions that allow for at least first-order 
effects of interrelations. 
A given system evaluation applies only to the particular situation being considered. 
Since the system may function well for some situations, only moderately for others, and Tl 
poorly for some, evaluations should be made for at least one situation of each of the types that ni 


are believed to be important. In this manner, an over-all evaluation of system performance 


can be obtained. co 
th 

THE MODEL d 
e 


For each time and basic component combination, the observed worth is expressed as 
the product of three quantities, which are so chosen that each factor has an intuitive interpre- 
tation. One of these quantities is a function of what happened and what should have happened 


for this basic component at this time. If the component performs accurately, this function has 
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the value unity. It has the value zero if the component does not perform at all. If the com- 
ponent performs so inaccurately that more harm occurs than if it had failed to perform, this 
function has a negative value. 

Multiplying the worth function of what happened and what should have happened are two 
nonnegative relative-worth constants. One of these factors is a measure of the relative worth 
of the component considered as compared with the other basic components for that time; the 
relative worths for the components associated with a fixed time are normalized so that their 
sum is unity. The other factor is a measure of the relative worth of the given time as com- 
pared to the values for the other times considered; these relative-worth values for the times 
considered are also normalized so that their sum is unity. 

Thus the measure of effectiveness, which is a sum of the observed worths, has the 
property that it is unity if the system is behaving ideally for the entire evaluation exercise 
and is zero if all of the basic components fail to perform throughout the exercise. Values 
near unity indicate good performance, while values near or below zero indicate poor perform- 
ance . 

Now let us consider the mathematical model for a specific evaluation exercise. Sup- 
pose there are I different times that are considered to be important. Associated with the jth 
time are Jj events (one for each component that is considered to be basic at time i), where 
what happens for each event is represented by a single number. 

The worth factor that depends on what happened and what should have happened for the 


jn event of the fh time is represented by the function 


F.= FAY... X,,): 
ij ha 5 has 


Here Vij is the number that was actually furnished, and Nij is the number that should have 


been furnished (by convention, Vij = -1 if no number is furnished). The factor Fij satis- 
fies the conditions 
Fij = Ly if Vij = Nij 
0, if V.. = =1 
1) 


<0, if Vij harmful. 


That is, this factor is unity if the accurate number is furnished, is zero if no number is fur- 


nished, and is negative if the number furnished is more harmful than no number at all. 


The worth factor that allows for the relative value of the j event of the jth time, as 


compared to the other Ji-1 events considered for this time, is denoted by Wii: The factor 


th 


that allows for the relative value of the i time, as compared with the other I-1 times, is 


denoted by Wi Both of these factors are nonnegative and normalized. That is, 


Ji I 


Wij 79, wi >0, >» Wij = 1, >» wil. 


j=1 i=1 
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Thus the observed worth for the jth event of the jth time is expressed in the form wjWijF ij; and 


the measure of effectiveness is 


The conditions imposed on the w,, Wij: and Fij imply that E=1 when Vij7Nij for all i, j; also 
that E=0 when Vij=-1 for all i, j; and that E can be negative if many harmful inaccuracies 


occur in the system operation. 


DETERMINATION OF INPUTS 

The inputs for the evaluation model consist of the times of importance, the basic com- 
ponents for these times, the values of the Vij the values of the Ni} and the values of the 
Fi j= F ij(Vij Nip» the Wij? and the w;. Great reliance on experience and judgment is needed 
in the determination of all these inputs. 

Ordinarily, the inputs are decided on after an examination of the recorded results of 
the evaluation exercise. Of course, specification of the times of importance, the basic com- 
ponents for these times, the functions Fij» and the factors Wij: w; prior to examination of the 
results of the exercise is desirable from one viewpoint; namely, subjective biases that can 
arise from such an examination are eliminated. However, this advantage seems to be small 
compared with the advantages of deciding on these items after examination of the exercise 
results; also the possibility of such biases can be reduced by the a priori adoption of some 
general nature guides for specifying these quantities. 

If every possible time of importance and the most detailed breakdown into basic com- 
ponents were used, evaluation of an exercise could represent a huge effort, even with the 
simplified model adopted. Experience with operations of similar systems combined with 
examination of the recorded results of the exercise should enable a large reduction to be 
made in the number of times of importance and in the number of basic components for a 
given time without much loss of sensitivity in the measure of effectiveness. 


One of the principal roles of the F ij? W.., and w; is to at least partially allow for 


important interrelations. Since an evaluation a of the type considered can usually take 
a huge number of fundamentally different courses, a priori specification of these quantities 
for all important possibilities does not seem to be feasible. Once the exercise is finished, the 
course followed can be determined and the Fij Wij: Wi specified for that situation. Also, the 
function F ij needs to be evaluated only for that Vij? Nij combination which actually occurred. 
The allowable values for the Nij do not include -1 but otherwise can be chosen in any 
suitable fashion. The actual values for the Nij and Vij are determined by examination of the 
results of the exercise. Here values for the Vij that are neither -1 nor any of the possible 
values for the corresponding Nij can occur because of gross errors in the system. 
Determination of the Fij leans heavily on experience, judgment, and the recorded 


results of the exercise. Only a single number is to be evaluated for each i, j combination 
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(i.e., the value for the Vjj and Njj values that actually occurred). However, in making this 
evaluation, taking the general shape of the function into consideration should be helpful. A 
graphical approach is frequently suitable for establishing the general shape of this function; 
sometimes the use of an elementary analytical function is satisfactory for this purpose (linear, 
exponential, etc.). The fact that Fjj is unity when Vjj=Njj and is zero when Vjj=-1 should be 
helpful in deciding on the general shape of the Fjj function. Special care is needed for cases 
where the observed value of Vjj is neither -1 nor any of the possible values for the corre- 
sponding Njj- 

Specification of the values for the Wjj and the wj also relies heavily on experience, 
judgment, and examination of the recorded results of the exercise. First- and second-order 
implications of inaccuracies should probably be examined for aid in the determination of 


values for these factors. 


DISCUSSION OF METHOD 

If a great amount of knowledge were available about a system under various operational 
conditions, a model might possibly be developed that would lead to a precise evaluation. Un- 
fortunately, this amount of knowledge is virtually never available for a new man-machine system 
of a complex nature and likely will not become available even after the system has been opera- 
tional for several years. Consequently, no "ideal" solution should be expected for this evalua- 
tion situation. Instead, the problem is to devise a feasible evaluation method that is capable of 
determining when a new system is functioning poorly and when it is functioning very well. Inter- 
mediate performance levels can be determined only in a rough sense because of experimental 
variation, lack of generality of the model, inappropriate choice of worth functions, etc. 

As already mentioned, the lack of knowledge about a system can be at least partially 
offset by careful use of such experience, judgment, and experimental type of information as 
may be available. Although somewhat subjective, information of this nature is the best that is 
to be had. While no direct experience with a system may be available, much data and experi- 
ence may be obtained about systems of a similar type. This information, combined with the 
judgment of personnel having appropriate backgrounds and the recorded results of the exer- 
cise, should be valuable in making the decisions necessary for the evaluation of the measure 
of effectiveness. 

The motivation for using a linear model falls under three headings. First, the avail- 
able information about the system performance ordinarily does not warrant a more sophisti- 
cated type of model. Second, a linear model is easy to comprehend; this is a valuable property 
in view of the strong part that judgment plays in deciding on the inputs. Third, a linear model 
is easy to handle from a computational viewpoint. 

The evaluation model presented is applicable to any specified subsystem as well as to 
the entire system. As a routine for application, the entire system should probably be evaluated 
first in a somewhat gross manner. If the measure of effectiveness has a value near unity, no 
further evaluation of this exercise should be necessary. If the system is not behaving 


Satisfactorily, more detailed evaluations can be made for the parts that are thought to be 






































18 I. M. GARFUNKEL AND J. E. WALSH 





performing poorly by using the results for the same exercise). This breakdown of evaluations 


is continued until the subdivisions which are causing trouble are isolated. 


PRODUCTION SCHEDULING EXAMPLE 

To illustrate how the evaluation model might be applied, an examination is made of 
the problem of assessing how well a manufacturing system is meeting a production schedule. 
Although not of a detailed nature or entirely realistic, this examination should furnish an 
indication of some of the considerations involved in adjusting the method to the situation being 
evaluated. 

A production process is considered to consist of separate basic operations, each of | 
which results in either a final product or in a product that is used in some further production 
operation. That is, a basic operation consists in taking one or more input products and manu- 
facturing an output product. Ordinarily, a given basic operation is carried out many times, 
furnishing repetitions of its output product. 

A production schedule is taken to be a specification, for each operation that is con- 
sidered to be basic, of deadline times for the production of the various repetitions of the 
output product of this operation. That is, for each basic operation, a production time limit is 
given for each of the outputs of this operation. A production schedule is usually stated in 
terms of some discrete time unit. For example, many schedules use a day as the time unit 
for specifying deadlines. 

To perform the evaluation of how well a manufacturing system is conforming to a new 
schedule, or meeting one that is already in operation, production is carried out with this 
schedule for some specified period of time. The product completion deadlines that fall within 
this time period represent the production schedule considered in the evaluation exercise. 

For the evaluation model, it is both convenient and meaningful to associate the times 
of importance with the deadline times that occur during the period of the exercise. Let the 
different deadline times be represented by the elements ti of a monotonically increasing 
sequence in which 1 < i <I. Then i=1 represents the first deadline specified by the schedule, 
and t is its time of occurrence; likewise i=I represents the final deadline specified by the 
schedule, and u is its time of occurrence. Here it is to be noted that the same deadline time 
can apply to several different product completions. For given i, the basic components are 
taken to be the output products which have the completion deadline t- If the discrete time unit 
used in stating deadlines is of a coarse nature, the value of Ji the number of basic compo- 
nents associated with the jth time, will usually exceed unity for large-scale systems. 

The value for Vij is taken to be the actual time of completion for the j™ product that 
has deadline time t. The value of Nii is t for all Ji of the basic components associated with 
the i*” time of importance. If some of the Vij are undetermined at the end of the time period 
specified for the exercise, the exercise is continued until these times are known. However, 
no additional deadline times or products are considered during this extension of the exercise 


period. By convention, time zero is taken as the time that the exercise begins. 
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In the case of production scheduling, obtaining the value of Fij = Fij (V ij? Nij) does 
not seem to be exceptionally difficult. Ordinarily, this function can be considered to depend 
only on the difference of its arguments (i.e., on Vij-Nij) and to have the value unity when 
Vij-Nij< 0. For Vjj-Nij > 0, the function Fij should nearly always be a monotonically decreas- 
ing function of Vij-Nij that can be graphically approximated on the basis of judgment applied 
to the available information. 

If a product whose deadline time falls within the exercise period happens to be com- 
pleted prior to the exercise period, the value of the Vij for this product is taken to be zero. 
Then the corresponding Fjj value is unity, as it should be for a product that is completed 
before its deadline. The case of Vjj = -1 (i-e., no number furnished) is not considered since 
a value should always be obtained for Vij: 

Trouble shooting is a standard component of nearly all manufacturing systems. Con- 
sequently, the extent of the difficulties arising from a given product failing to meet its dead- 
line, as compared with the difficulties arising for the other products of the system, should be 
at least roughly known. This general knowledge, combined with the actual results of the 
exercise, should furnish a reasonable basis for specifying the values of the Fij; Wij; and wj. 

The values of I, the Jj, the Fjj, the Wij; and the wj, being given, evaluation of the 
measure of effectiveness is straightforward. If E is near unity, the system is conforming 
very closely to the given production schedule. If E is near zero (or negative), the system is 
falling far behind the specified production schedule. Intermediate values of E represent inter- 
mediate levels of conformance to the production schedule evaluated. The quantitative impli- 
cations of intermediate values of E depend on many considerations, such as the forms of the 
Fij functions, the values for Wij and w;, and the appropriateness of these quantities for rep- 


resenting interrelations. 
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This paper essentially describes the activities of the RAND Logistics | 


Systems Laboratory during its first year of operation. During this year, 


the laboratory conducted an important first experiment in which it tested | 


a series of proposed changes to the Air Force supply system. 


The test involved the creation of two versions of an Air Force logistics 
system, one which simulated a 1956 Air Force logistics system and an- 


other whichwas modified to contain the proposed supply policies. These | 


two systems were operated in parallel and were given the same support | 


task to accomplish. The system which was able to accomplish its job at 
lesser cost was deemed the better system. 


The simulation was a composite of man-machine models, in which the 


men were assigned management responsibilities and functions and the | 


machines performed routine actions. The men were identified by titles 
which have counterparts in the real-world logistics system. The two 
logistics systems were operated under compressed time so that each 


hour in the laboratory simulated a day's experience in the real world. | 


This enabled the simulation to cover 3-1/2 years' experience and 2 | 


simulated wars. 


The experiment indicated that, according to the criterion established, | 
the logistics system containing the newer supply policies was a better | 


system. The experiment also indicated ways in which revisions to the 
procedures by which the new policies are implemented would benefit 
their future use. 


The experiment indicated that this type of logistics system simulation 
did develop many more insights into the operation of proposed policies 


than traditional analytical techniques normally provide, and it did extend | 
the traditional descriptions provided to real-world organizations of such | 


proposed policies. At the same time, the experiment did highlight the 
complexities of testing logistics policies, even under the idealized con- 


ditions of the laboratory, and the problems of interpreting experimental | 


findings. 





INTRODUCTION 


the Air Force with additional information and data about logistics research findings than are 


obtained through traditional research techniques. 





+ 
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The RAND Logistics Systems Laboratory was begun in October 1956 to help provide 


The basic technique to be used in the 
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Laboratory was to be man-machine simulation experiments. Although such simulation had 


been used previously, it was designed for different purposes, therefore, the Laboratory was 





essentially required to develop a new method. 

As a means of obtaining experience with logistics simulation by use of man-machine 
techniques, the Laboratory started by running two exercises early in its life!» * (hereafter 
referred to as Prologs 1 and 2). These served as a prelude to Laboratory Problem 1 (LP-I), 
which is the subject of this paper. 


EXPERIMENTAL DESIGN 
Objectives of LP-I 

One objective of LP-I was to examine, for feasibility, a series of proposed changes to 
the logistics system. The modified logistics system containing these changes was called LS-2 
(Logistics System 2) during the experiment. Another objective was to compare the perform- 
ance of the proposed policies with an alternative but an aribtrary system, called LS-1 
(Logistics System 1), which was intended to simulate the 1956 Logistics System of the Air 
Force. The design used in the experiment for comparing the two logistics systems was to 
have both logistics systems perform the same job in terms of numbers of bases, aircraft, fly- 
ing hours, etc., supported and then to measure the differences in cost and aircraft readiness 
resulting between the two systems from their use of different logistics policies. Under this 
design, the better system is that one which does the specified job at lesser system cost or 
greater aircraft readiness or both. 3 


Realistic Simulation 

To test adequately the proposed policies, it was necessary to reflect in the simultion 
the complexity of the logistics system with which the policies would have to deal. It was also 
desired to provide a realistic setting for the experimental participants who would use these 
policies, so that their experience would simulate what might occur in the real world. In order 
to reflect the complexity of the logistics system, and to provide this realistic setting for the 
participants, it was necessary to include considerable detail of the logistics environment in 


LP-I, as the succeeding discussion will indicate. 


Total Weapon System Life 
The proposed changes to the logistics system used in LS-2 were intended to function 


over an entire weapon system life and to produce savings in system cost over this lifetime. 


1 Logistics Systems Laboratory, ''First Tooling-Up Exercise for Logistics Systems Laboratory 
(Oct. — Nov. 1956)", The RAND Corporation, RM-1924, July 1, 1957. 

- Logistics Systems Laboratory, ''Second Tooling-Up Exercise for Logistics Systems Labora- 
tory (Jan. — Feb. 1957)," The RAND Corporation, RM-1961, August 19, 1957. 

3 The official RAND report is R-323, ''Laboratory Evaluation of Supply and Procurement 
Policies: The First Experiment of the Logistics Systems Laboratory.'' As a reading of this 
paper will indicate, this type of simulation is a large-scale group activity, and it is the 

quality and effort of the group that determines the outcome of the experiment. 
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Consequently, one of the requirements of the experiment was to study the behavior of the pro- 
posed policies over the total weapon system life. 

These proposed policies had to be able to deal with the conditions and problems that 
exist at the phase-in of a new weapon as well as those that occurred at the phase-out of the 
weapon. Such conditions in the real world are oftentimes dynamic and uncertain in behavior, 
because they include changes in aircraft series and engineering modifications and revisions in 
Air Force programs as defined by numbers of aircraft, flying hours, etc. In addition, the 
policies should be able to utilize past history as the weapon system experience unfolds. 

The weapon system life for a fighter-type aircraft, which was simulated in LP-I, lasts 
for about five years. To attain so much experience in the three months planned for the running 
of the experiment, it was necessary to effect significant time compression, so that, for 
example, a week of real time would simulate several months of simulated time. Time com- 
pression was further produced by revising the time units so that a laboratory month contained 
only ten days instead of the thirty days of the real world. Such a time-unit distortion was also 
used in Prologs 1 and 2, and analysis of their results indicated that this change did not 
seriously bias or damage the data produced by the simulation. 


Two Logistics Systems in Parallel 

To attain the objective of comparison between two logistics systems, it was important 
that the two systems function in parallel as simulated time transpired in the laboratory. This 
was done to facilitate experimental control of the two logistics systems by insuring that they 
were both treated the same. For example, policy questions or interpretations were frequently 
required during the experiment, and it was important that both systems be equally affected by 
these interpretations. By knowing how each system stood at the same point in time, it was 
much easier for the experimental staff to make decisions that would not prejudice one or the 
other systems. Also, these decisions could be introduced with a minimum of disruption to both 
systems since they did not have to be ad hoc for one system and post hoc for the other. 

In operating both LS-1 and LS-2 in parallel, we viewed LS-1 as serving, first, to pro- 
vide a bench mark for realism since it was intended to simulate the 1956 Logistics System of 
the Air Force and, second, to provide a standard against which any benefits realized with the 
proposed system could be measured. It should be emphasized that LS-1 is not the present Air 
Force system, and, as a matter of fact, in view of the complexity of the Air Force logistics 
System and the different ways it operates in the hundreds of organizations comprising the 
logistics system, it would be very difficult if not impossible to state that any single simulation 
represents the Air Force logistics system of any date. Additionally, we have the further con- 
dition that the policies used in LS-1 (the 1956 Air Force Logistics System) differ from the 
policies used in the Air Force today, which is to be expected in a dynamic organization. 

LS-2 contained limited modifications in the logistics system, particularly in supply 
policies. These proposed changes were the result of RAND and Air Force research over the 
past several years. Some of these proposals are now being introduced into the Air Force. The 
experiment was intended to provide an opportunity to examine how the complete set of policy 


changes would function in a revised system if they were fully implemented. 
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Aircraft Parts Sample 





The simulated aircraft in LP-I was an ADC fighter interceptor of current vintage. In 
order for the volume of data created in the simulation to be handled feasibly, the aircraft was 
represented by a sample of 800 parts drawn from its entire fuselage. As a further simplifi- 
cation, this sample included only peculiar airframe parts. Consequently, the experiment did 
not provide the opportunity to study the interaction among weapon systems in their demand for 
common spare parts, which is not a trivial problem to the real-world logistics system. 

The specific sample of parts used in LP-I was biased towards higher demands to 
exercise the responsive character of the policies. Furthermore, the sample was picked so 
that it contained a higher proportion of hi-valu (higher priced) parts and a lower proportion of 
cheap parts. This was done to provide enough sufficiently accurate data to evaluate the per- 
formance and cost of the hi-valu policies. However, cheaper parts predominated numerically ( 
in the sample, so the effects of such parts on logistics system cost and performance were still 


importantly represented. The extent of bias in the sample as to the demand rates and prices of ( 


its parts was known, so it was possible to adjust the results for such biases, as will be seen later ‘ 
t 

Failure Model i 
An important triggering action in a logistics system that supports aircraft is the fail- t 

ures which occur on the aircraft as a result of flight, age, inspection, etc. The validity of the I 
results in the experiment depend significantly on the realism with which aircraft failures canbe V 
simulated in the laboratory. The mechanism for generating failures during the simulation was I 
called the failure model.* For this reason, much real-world data which provided information i 
on the pattern with which failures occur within particular assemblies, at particular Air Force base 
locations, and over the different phases of the weapon life were represented in the failure model. € 
i 

Man and Machine Interaction i 


In LP-I considerable use was made of electronic computers and other mechanical com- 


puting devices. This was necessary to handle the volume of data and the time compression i 
used in the experiment. An important part of the design was to decide which of the simulated k 
logistics operations were to be done by machines and what decisions and actions were to be ¥ 


assigned to the human beings participating in the experiment. 


The guiding rule finally chosen was to use the computers to simulate the physical actions 
that occur within the logistics system, such as the inspection of aircraft, the removal of de- , 
fective parts, the replacement with serviceable parts, etc.; the thousands of clerical actions . 
which are necessary to maintain records, prepare reports, fill out forms, etc.; and the F 
thousands of routine decisions which are made by lower-level management personnel who follow 
closely prescribed rules and procedures. Thus, the actions left for the human beings in the R 


experiment were those which involved major management decisions for which routine rules had 





unl 


4Sce W. McGlothlin, G. C. Noonan, Jr., G. E. O'Dell, and S. L. Pollack, ''The Simulated 
Aircraft and Its Failure Model in LP-I,'' The Rand Corporation, RM-2177, May 26, 1958. 
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not been prescribed or which affected the operation of the logistics system in a very significant 
way. Necessary reports and data to be given the human beings for making such important de- 
cisions were a basic part of the procedural design which went into LS-2, and the formats of 
the regularly prepared reports of the Air Force were used for LS-1. Thus, the reports and 
data used for decision making in LS-2 may have much relevance for implementation of that 
system in the real world. 

This decision to confine the human being to making only important decisions was partly 
determined by the experience in Prolog 2, which pointed strongly to the boredom and loss of 
morale which could occur in such experiments if the human subjects were used to make re- 


petitive and straightforward decisions in experiments lasting several weeks. 


Organizations 

The Air Force organizations simulated in LP-I included: Headquarters USAF, Head- 
quarters ADC, Headquarters AMC, the weapons system manager's office, several air bases, 
anda factory. These organizations can be classified into two groups: the embedded organiza- 
tions, which included the higher headquarters and the factory and the floor organizations, which 
included the bases and the weapon system manager's office. The policies under study were 
those which are made primarily in the weapon system manager's office, since they deal with 
procurement and distribution of spare parts, as will be discussed soon. For that reason, the 
weapon system manager's office was the central organization in the study. The bases were 
peripheral organizations which were simulated in sufficient detail to provide reasonably good 
inputs to the weapon system manager for his decision making and actions. The embedded 
organizations were manned by laboratory staff personnel because their primary function in the 
experiment was to provide a surrounding environment for the floor organizations and to produce 
in this environment the kinds of situations and conditions which tested both logistics systems 
identically. 

The bases were assigned real names, and they were assumed to be located at these real 
places so that the effects of transportation delays and differences in transportation cost would 
be part of the environment in which the logistics policies had to function, as is true in the real 


world. 


POLICIES 

The policies tested in LP-I were: responsive procurement of hi-valu (high-cost) items; 
economical procurement, repair, and distribution of (medium and low-cost) Category 2 and 3 
parts;5 and automatic resupply and centralized record keeping. These policies will next be 


discussed briefly in turn. 


Responsive Procurement of Hi-Valu Items 
Historically, the Air Force has provisioned for many airplanes before the first opera- 


tional aircraft has been delivered from production. At this early time in the life of a new 


eee 





5 

Category 2 parts are those with unit prices between $10.00 and $500.00, not classified as 
hi-value parts for reasons of high demand. Category 3 parts are those with unit prices less 
than $10.00. 
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aircraft, there is great uncertainty as to future requirements for spare parts, both because of 
difficulties in predicting failure or demand rates and because of anticipating program changes. 
An approach to reducing the amount of uncertainty with which the logistics system has to con- 
tend in procuring spare parts® might be to use an institutional solution in which the provision- 
ing of hi-valu parts would be deferred as long as requirements of the system permit. 

Under deferred provisioning of hi-valu parts, during the phase-in period, the logistics 
system would depend primarily on a buffer stock established at the factory, from which it 
would draw as demands occur within the Air Force. By living primarily out of this buffer 
stock during the phase-in, it might be possible to defer the final procurement decisions on hi- 
valu parts by one to two years. In this way, the Air Force might be able to live through a very 
uncertain period in the program of a new weapon system without unnecessary procurement of 
hi-valu parts. 

The importance of improving provisioning for hi-valu items is indicated by the fact that 
40 to 50 percent of the money spent on aircraft spares is spent on hi-valu items, although the 
hi-valu items number less than 5 percent of the total line items for an aircraft. 

This policy is attractive because the demands for spare parts are typically very low, 
and therefore a small buffer stock at the factory provides considerable protection for most of 
the spare parts. For those few parts with high demands, the buffer policy would be replaced 
by a formal procurement as soon as adequate experience were available to estimate the demand 
rate accurately. For the low-demand items, deferment could be practiced until the final de- 
cision on procurement had to be made before the production line shut down. In most cases the 
buffer would be more than adequate for these line items, and so no further procurement would 
be necessary. It might even be possible to use a portion of the buffer in the final aircraft 


assembly.’ 


Economic Procurement, Distribution and Repair of Category 2 and 3 Parts 

Until recently, the cost of ordering cheap parts was higher than the value of the stock 
being ordered. This was especially true of Category 3 items. Despite their low price, these 
items affect importantly the performance of the logistics system, because historical studies 
have shown that Category 3 items account for about 40 percent of aircraft out of commission 
for parts (AOCPs). 

A more comprehensive formula has been developed at RAND for computing the procure- 
ment, distribution, and repair requirements of Category 2 and 3 parts.® This new formula 
takes account of many more factors than those considered in current Air Force policies on 
procurement, distribution, and repair of Category 2 and 3 parts. For example, in addition to 


the price and demand rate of the part, the formula considers as additional factors the cost of 


This is especially important for the hi-value spares. 

J. Petersen, 'Savings from Procurement Deferral with Interim Contractor Support; The 
Case of High-Value Airframe Spares,'' The RAND Corporation, RM-2085, January 10, 1958. 

F. Ferguson and L. Fisher, ''Stockage Policies for Medium and Low-Cost Parts,'' The 
RAND Corporation, RM-1962, April 18, 1958. 
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reordering, the cost of holding inventories, and the cost of parts shortages. The effect of 

this new formula is higher reorder points and reorder quantities for the cheaper parts than 
provided by previous policies. These higher levels increase the protection provided by supply 
inventories at a base, thereby reducing the number of aircraft out of commission for parts 

and the number of priority requisitions required. The lesser requirement for priority action 
reduces the cost of processing requisitions at base and depot and of transporting cheap parts 
to the base. The intended net effect is therefore fewer shortages at practically no change in 
cost. This, then, was the second new supply policy built into LS-2. 


Automatic Resupply and Centralized Record Keeping 

In the past, bases have computed their own stock control level and done their own re- 
ordering. The effort involved in this inventory control was directly proportional to the number 
of line items stocked. Therefore, the tendency was to limit the number of line items stocked 
at the base.” However, no explicit systems were given for determining how to select those 
items to be stocked. This situation accounted somewhat for the comparatively high proportion 
of AOCP's caused by cheap items. 

Also, under the present Air Force system for inventory control of Category 3 items, 
once such stock is shipped from the depot it is assumed to have been consumed and therefore 
no longer subject to inventory control. This policy leads managers to keep more stock of 
such items at the depot so that the stock can continue under their control. This practice con- 
flicts with the apparently desirable policy of stocking bases heavily with cheap parts. 

The introduction of electronic data-processing equipment might provide an opportunity 
to alter this pattern, especially for Category 3 items. Under the concept of the data process- 
ing center (DPc),!° this center would maintain up-to-date knowledge of stocks in all cost 
categories by location. The center would adjust this knowledge through frequent reports 
received from the using activities so that current inventory information would always be 
available. The DPC would also compute the stock control levels and reorder points, and 
thereby it would effect automatic resupply. Thus, by improving management control of the 
stock, such DPC's might help greatly to reduce the number of AOCP's suffered in the Air 
Force as well as the total Air Force investment in stock. LS-2 contained a DPC which 
performed centralized record keeping, recomputation of stock control levels and reorder 
points, and automatic resupply. 

Although the DPC simulated in LP-I did provide much information about the operation 
of such a facility, it was not a complete simulation of a real-world DPC. For example, the 
Supply transactions reported to the DPC were completely generated within the computer. 
Therefore, the clerical errors in reporting transactions that would occur in the real world 
were not encountered by the simulated DPC. Also, all transactions were reported individually 





9 The Strategic Air Command initiated SACSIP (SAC Supply Improvement Program) in 1954 
and thereafter set ceilings on the range of line items each type of base could stock. Other 
commands followed it. 

10M. A. Geisler, ''Some Principles for a Data-Processing System in Logistics, '' Naval Re- 
Search Logistics Quarterly, June 1958. 
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to the DPC, regardless of price or any other characteristic of the line items involved. This 
may not be feasible in the real world because of the volume of reporting that would be involved 
with individual transactions. For cheaper parts, a summary form of reporting might be used. 
Furthermore, there was no distinction in the simulation between the paper records that would 
be kept at the DPC and the physical inventory that would be located throughout the logistics 
system. Both were represented by the same set of numbers in the computer. Therefore, 
discrepancies between records and inventory did not occur in the simulation, although they 
would undoubtedly occur in the real world. These simplifications were introduced to reduce 
the amount of machine programming and the volume of machine calculation which would other- 

wise have been necessary. At the same time, these simplifications were not believed to have 


reduced the validity of the basic policy comparisons between the two logistics systems being 


studied. 

In summary, the three sets of policies just described were what characterized LS-2. I 
In LP-I, we wanted to test the feasibility of using these policies in a single logistics system V 
(LS-2) and compare their performance and cost with the policies used in LS-1. a 

[ 

PREPARATIONS FOR THE LABORATORY RUN t! 

An indication of the amount of work required to prepare for LP-I may be gained from t 
the following statistics. First, the machine programs, which were used to simulate many of 0 
the elements represented in the logistics system, included about 25,000 different machine tl 
instructions. Furthermore, over 30 participants were used on the laboratory floor to act as si 
base managers and as personnel in the weapon system manager's office in LS-1 and LS-2. In 
addition, about 50 persons were on the laboratory staff. In order to train the participants for hi 
the on-floor activities, two manuals (one for LS-1 and another for LS-2), each containing it 
about 350 pages, were prepared for the experiment. Also, approximately 100 forms were le 
used during LP-1 to simulate the different forms and reports required to manage the simu- as 
lated logistics systems. Some of these forms were standard Air Force forms, but many of pe 
them were new forms prepared to fit the new policies and the special needs of the laboratory m 
activity. Finally, each of the participants was given from 2 to 4 weeks training so that he 
would be familiar with the policies and procedures which he was expected to use during the 
experiment. 
OPERATIONS OF LP-I " 

The main elements of the simulation consisted of the provisioning conferences held : 


before any bases were activated and the daily routines and quarterly activities of the base and 
cle 
weapons system manager. Each of these situations contain important features of a real- 


world logistics system. In the remainder of this section the procedures used in each of these 
situations will be described, as well as the supporting laboratory activities, including the oe 


1] 
organization of the statistical operations and analysis. 
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The Provisioning Conference 





i d As is true with the initiation of a new weapon system in real life, LP-I began witha 

- provisioning conference for LS-1 and LS-2. The participants at the provisioning conference 

ald made decisions very much like those made in the real world, that is, they established pro- 
curement and delivery schedules for spare parts, and they determined the initial tables 
representing amounts to be stocked at the first few bases at which the new weapon system 
would be assigned.! : 

ye The Bases 

ner- LS-1 Bases 

ave The activities of the LS-1 bases were primarily simulated by means of a computer 

ng program which calculated the daily experience of the base as follows: 

Each day, the aircraft at the LS-1 bases tried to fly the flight schedules that were 

2. provided to them by the experimenters.!2 As a result of this flight activity, malfunctions 

2m were detected which required maintenance personnel and spare parts. The repairs on these 
aircraft were accomplished as rapidly as possible within the resources available to the base. 
Daily reports were given to the base manager, an assigned Air Force airman, which told him 
the aircraft status, flights accomplished, and spare parts needing special management atten- 

“om tion. The management decisions were then communicated through the computer, or by mail 

y of or telephone, to the appropriate organization for action. The base manager's objective in all 

2 this was to satisfy his assigned alert and flight program, and he had to obtain as many re- 

Las sources as necessary to do so. 

In At the end of each quarter, the base manager received a balance listing which gave 

s for him the current inventory of each line item at his base and the number of units of each line 

4 item issued during the past quarter. He used this information to recompute his stock control 

“e levels and reorder points and to determine the changes in his inventory which should be made 

1u- as a result of these recalculations. He also reviewed his requirements for field maintenance 

r of personnel for the next quarterly period and determined whether or not a change in his field 

tory maintenance allowance should be made in view of his anticipated field maintenance workload. 

ws LS-2 Bases 

on The LS-2 base managers were more limited in the management activity they could 
perform because their inventory levels and records were maintained centrally by the Data 
Processing Center, and the weapon system manager's office established the inventory levels. 

id However, these managers did assign priorities to those parts that were repaired at their 

iil base. During the simulation, the lessened responsibility of the LS-2 base managers was 
Clearly evident. 

these 

he Sy 


1] aes : ; : : , : ; 
Two series of aircraft were simulated in the experiment, with the D series differing from 
the C series in some majo r aspects. 
Lp ; ‘ : ; ; P 

in effect, the flight program was exogenously introduced. This was done to help insure 


that both systems met the same flying program. 
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next quarter. 


The Weapon System Managers 
The LS-1 Weapon System Manager 





aircraft sent to him periodically from the bases. 


from the computer reports of work accomplished and shortage problems. 


sulted in corrective actions which were then communicated to the relevant 
this reality. 


and consumption report similar to that used in the real-world Air Force. 
the inventory balances, issues, and due-ins of the Category 1 and 2 parts. 


parts repair manager should follow for its next quarter. 


13 Inspection and repair as necessary, performed at depot. 


seems appropriate for very large numbers of low-cost items. 


At the end of each quarter, the LS-2 base managers received balance listings like 


two constituted his supply organization and the latter two his maintenance organization. 


those given to the LS-1 managers. They used these listings primarily to follow the operation 
of their bases, since they had no responsibility for the computation of inventory levels or 
routine requisitioning. These latter actions were performed by the Data Processing Center. 
They also used this information to estimate what their field maintenance workloads would be 


and to determine if they had the requisite field maintenance personnel at their bases for the 


The LS-1 weapon system manager was assisted by a requirements manager, a dis- 
tribution manager, an IRAN?? manager, and a parts repair depot (PRD) manager. The first 
These 
5 people were experienced Air Force personnel. Their function was to keep the bases supplied 
with parts, accomplishing this task at minimum cost. In addition, the IRAN manager repaired 


Here, too, parts of the organization were simulated on the computer, including the 
maintenance and repair shops. Daily, the weapon system manager's organization received 
These reports re- 
organizations. The 
intent in all this activity was to simulate in appropriate detail the kinds of logistics problems 


that such real-world managers would encounter and therefore to have their decisions reflect 


At the end of each quarter, the LS-1 requirements manager received a stock balance 


This report covered 


For Category 3 


information, the requirements manager relied on the balance listing of the depot of Category 3 
parts.!4 He also received a Supportability Cost Analysis, which is a quarterly requirements 

calculation specifying the procurement and repair schedules. This report was based upon the 
procedures used at the real-world Sacramento Air Materiel Area. The quarterly support- 
ability calculation gave the requirements manager all changes in procurement and master 

repair schedules required since the previous quarterly report. Using such quarterly reports, 
the requirements manager was able to determine what orders should be placed with the factory. 


He also gave to the parts repair depot manager the changes in the schedule of repairs that the 


14The assumption behind using this report is that depot issues of Category 3 parts are accu- 
rate enough measures of consumption. Although there are shortcomings in using sucha 

report, for instance, depot issues also include requests by bases for inventory to support 
stock levels, the use of depot issues is a much simpler and cheaper reporting technique which 













































ite 
sc 
an 
se 


Sir 






on 


lis- 
first 
these 
plied 


aired 


The 
ems 


lect 


ince 
overed 
y3 
gory 3 
nents 
on the 
rt- 
aster 
sports, 
factory: 
hat the 


accu- 
a 

ort 

e which 





A FIRST EXPERIMENT IN LOGISTICS SYSTEM SIMULATION 31 


The LS-1 distribution manager used the Stock Balance and Consumption report and the 
palance listing of the depot for Category 3 items to recompute the various inventory control 
levels required for the management of his storage site. The distribution manager also deter- 
mined whether or not redistribution of stocks, especially of Category 1 and 2 parts, was 
required among the bases of his weapon system. 

The IRAN manager received a quarterly report of the amount consumed of each line 
item used in IRAN. He used this information with his schedule of future IRAN's to determine 
the amount of inventory that should be at his activity during the next quarter. 

The PRD manager received a quarterly report of materials consumed and man-hours 
used in the repair of assemblies. He used this experience in combination with earlier experi- 
ence to determine his man-power and material standards for the next quarter. The PRD 
manager used these adjusted standards along with the master repair schedule to establish the 


PRD material and man-power requirements for the next quarter. 


The LS-2 Weapon System Manager 





The responsibility of the weapon system manager was organized somewhat differently 
from that of the LS-1 weapon system manager. He had assisting him a hi-valu manager, a 
Category 2 and 3 manager, an IRAN manager, and a parts repair depot manager. The hi- 
valu manager was in charge of all management activities concerned with the high-cost parts, 
while the Category 2 and 3 manager was concerned with the medium and low-cost parts. In 
LS-1 the responsibilities were split by function (requirements and distribution) and not by 
groups of parts. The IRAN and parts repair depot managers performed essentially the same 
functions in LS-2 as in LS-1, except that the inventory control and the routine resupply of 
their spare parts was handled through the Data Processing Center. 

Although the personnel filling these positions were experienced with Air Force systems, 
considerable training was required to teach them the new policies and procedures that 
characterized their logistics system. To an extent, they probably never learned the system 
completely, and RAND personnel had to be available for technical consultation. 

The hi-valu manager received a monthly report which indicated to him those hi-valu 
items that required procurement and repair. In using this monthly report, he followed pre- 
scribed mathematical rules to determine the orders to be placed with the factory and the 
amount to schedule with the parts repair depot manager. The report was monthly during the 
sensitive phase-in period, but as this period passed the report becamse quarterly in frequency, 


since the system did not have to be as responsive in the later phases of the weapon system life. 


The Category 2 and 3 manager received no quarterly reports; instead he reliedon daily 
data. Each day, he was told which parts required procurement and which parts required re- 


pair action, and he was expected to take this action at that date rather than waiting for the end 


FS 


of the quarterly periods. The technique also updated consumption rates, etc., as experi- 


ence accumulated. 





15 
F, Ferguson, and L. Fisher, "Stockage Policies for Medium and Low-Cost Parts,'' The 
RAND Corporation, RM-1962, April 18, 1958. 
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The Data Processing Center reported organizationally to the weapon system manager. 





































He prescribed the policies to be followed by the Data Processing Center in effecting require- 
ments and distribution actions. The bases were required to submit transaction data to the 
Data Processing Center. Obviously, this latter requirement would have to be met with the 
approval of Air Defense Command or Headquarters, USAF — the two present command head- 
quarters able to enforce relationships between the weapon system manager (AMC) and the 
bases (ADC). 


The Factory 
The factory was an embedding activity, that is, its actions were controlled by detailed 
rules specified by the experimenters. It was thus part of the environment in which the policies 
of both logistics systems had to function. a 
The factory manager issued quarterly procurement status reports to the weapon system 
managers of each logistics system. This was done primarily to provide the managers with 
accurate information of what was due in from production, since such data were critical for e 


accurate production decisions. 


The EAM and 704 Computer Activity 


Each simulated day in the experiment took one hour of laboratory time. Half of the I 
hour was required for computer operation, and the other half hour was used for the manager's : 
review of the computed reports, his determination of data to go to the computer, and the f 
preparation of the data for the computer. Thus, the simulation activity was completely n 
dependent on the computer operation.!® Therefore, whenever difficulties transpired at the ; 
computer because of either input-data errors, programming stops, or computer breakdowns, s «S 
the simulation activity ground to a halt. Because of these difficulties, the actual computer pt 
time required to complete a successful run probably was, on the average, close to an hour. th 


Therefore, since we ran approximately 400 simulated days in the experiment, about 400 


computer hours were required. In addition, each quarter's reports took an additional 4 com- - 5 
puter hours, so that altogether these reports required a total of 56 hours. It wouldnot bean | "| 
understatement to say that the total experiment required at least 500 hours of computation ina ; @ 
large-scale elect ronic computer. as ° 

ar 
The Control Point sk 


The control point handled all the data received from the floor personnel and prepared 


them for the computer. It also distributed the listings received from the computer. The con- T 
trol point also performed the highly essential function of checking all inputs to the computer 

for procedural consistency because input errors, especially in the early stages of the experi- of 

ment, caused many computer reruns. The control point also helped to simulate the mail = 

Th 

—— frc 
6A very large-scale computer was used that had several magnetic tape control units and a 

32,000 word core memory. i 





. manned by laboratory personnel. It performed a vital training and orientation function because 
it received many calls from personnel on the laboratory floor requesting explanations about 
the origin of numbers on various reports. Statistical Services also did much trouble shooting 
during the operation because many of the questions asked by the floor personnel required 

iled examination of the reports and investigation of the procedures. 

‘cies The day-by-day experience of Statistical Services in answering questions of the floor 
personnel and in repairing the data system when it was not functioning properly produced con- 

stem siderable system modification and improvement during the run. Thus, in many ways it was 

h the nerve center and the eyes and ears of the laboratory staff as far as the machine operation 

r and the accuracy of the data inputs and outputs of the simulation were concerned. 

The Observation Program 

An observation staff used a variety of systems during the experiment to gather be- 

e havioral information from the floor. A large and complex telephone and microphone recording 

ger's system was used to record as much as possible of the information and conversation that was 
passed orally among the floor organizations. Also, personnel from the observation staff man- 
ned stations on the observation deck, which permitted them to observe visually and aurally what 
he was happening on the floor and to make notes so that later studies of the recordings would be 
wns, supplemented by additional information concerning the situation that prevailed at any particular 
or time. These data were intended as an aid in understanding the interorganizational activity and 
ur. the conflicts arising between organizations as a result of the experiment. 

In addition, at frequent intervals questionnaires were administered by the observation 
com- staff to the floor personnel to determine whether or not the personnel felt the simulation was 
be an realistic, what problems the personnel were encountering in the simulated system, what 
on in a experiences they were having in using the policies, and their general satisfaction with the 

experiment. These questionnaires were supplemented by interviews held between observation 

and laboratory personnel in an attempt to determine the comprehension that the floor personnel 

Showed in their use and discussion of the policies and the procedures representing their system. 
pared 

— The Analysis Program 

suter An analysis staff also was used, which was concerned primarily with the performance 

xperi- of the policies in the experiment. The analysis staff maintained daily and current information 

il on the simulated logistics systems, especially on their comparative effectiveness and cost. 
These data were derived primarily from the daily and quarterly reports which were required 

. from the personnel on the floor. The continuing knowledge of the current status of each 
ind a 
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delivery, transportation, and repair time lags by holding the punch cards, which simwated the 


parts, for the requisite number of days required for the lag-time simulation. 


Statistical Services 


This organization, one of the most important in the operation of the experiment, was 


logistics system by the analysis group also helped much in quality control, because sharp 
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differences between the two systems were investigated and explanations sought. This activity 
was most invaluable in helping the laboratory staff to control the experiment and to assure that 
invalid conditions did not continue. 

The logistics system simulation included the monetary inventory accounting system 
used in the Air Force for maintaining dollar information on the inventory status for aircraft 
spare parts. This information was provided to the floor participants in a form very much 
like that used in the real world. The floor participants did not make much use of this informa- 
tion (as is the case in the real world), but the analysis staff found these aggregative data very 
useful in understanding what was happening to each logistics system and in determining the 
extent to which the floor participants were following the policies. 

The fact that the two logistics systems were being run simultaneously and in parrallel 
was a great help to the analysis staff in detecting discrepancies or other factors which tended 
to affect one system or the other unfairly or ina different way. By having the same data for 
the two systems at the same point in simulated time, we were able to establish quickly where 
the systems differed and then to determine whether these differences were those due to the 
policies or those due to the artificiality of the simulation, including possibly the individual 
differences of the managers themselves. 

Each quarter an analysis and observation report was issued which was given only to 
the laboratory staff to enable them to stay abreast with the results of the experiment and to 
suggest improvements, changes, etc., that were necessary to keep the experiment operating 


satisfactorily. These reports were an invaluable documentation of experience during the run. 


PEACETIME RESULTS 
Qualification of Results 

Before the results obtained in the experiment are presented, it is important that the 
circumstances surrounding the development of these results be understood. These qualifica- 
tions are developed in detail to point up the complexity of logistics experimentation and the 
complex nature of the phenomenon being studied. Despite these qualifications, the methods de- 
scribed in this paper seem to offer distinct benefits in logistics systems evaluations and study. 


Complexity of System 





Complex experiments of this type inevitably involve many difficulties. There are no 
apparently simple ways to attack these types of problems. Simpler simulations can be de- 
vised, but they would answer a different set of questions from those which this experiment was 
intended to do. 

The first factor to be considered is the complexity of the system. As has been indi- 
cated earlier, we endeavored to maintain realism in the simulation by introducing as many of 
the relevant elements as we could identify in the real world. With such a complex system, 
discrepancies that occur in one part of the logistics system and that may be due to misunder- 
standings of the simulation may spread throughout the system, and their effects on the data 
and results may be reflectedeverywhere. This means that each important finding in the experi- 


ment must be assessed as to whether they are true effects of the logistics system or artifacts. 
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vity With such a complex system, it was also very difficult to maintain the degree of 
that experimental control necessary to insure complete validity of results. For example, this 
included the objectives of having the aircraft in both systems fly the same number of hours 
' each day so that both systems performed exactly the same job. This was impossible to ob- 
ft tain, and some of the ad hoc changes in flight programs that had to be made to attain this 
objective may have affected the results somewhat, but not importantly. 
ma- It was also very difficult for the laboratory staff to be sure that the floor personnel 
pery were always following the prescribed policies. Whenever and wherever we could, we did try 
e to make sure that they did follow the policies; and by appropriate directives from higher head- 
quarters and by other means, the necessity to observe policy was constantly brought to the 
lel attention of the floor personnel. 
ided 
lie Procedural Difficulties 
Many procedural difficulties occurred during LP-I. One of the most important was the 
ais programming errors which had not been fully eliminated before the run began. Although a 
: little checking of the program was done before the start ofthe run, not enough time was allowed 
in our schedule for this activity. Very likely, it would have been impossible to have checked 
the programs completely because of their complexity, with their thousands of machine instruc- 
7 tions and involved logic. Thus some programming errors were present, and it is undoubtedly 
_— true that their effects are reflected somewhat in the results. 
Input errors made by the managers in submitting data for the computer also were de- 
ae tected; some of them we found before they affected the results, but for others we were too 
late. These errors also probably influenced the results obtained. 
Data Difficulties 
the Data difficulties were among the most important of the effects which color the results 
fica- of LP-I. These difficulties were largely in the failure (or demand) rates used for the spare 
the parts in the sample. A most important error was the bias in the hi-valu demand rates which 
ie ities was not discovered until the experiment was well under way. This happened because we were 
ale. given estimates of future consumption developed for use in computing future Air Force re- 
quirements rather than data on past Air Force experience. The effect was that the demand 
rates for hi-valu items used in LP-I were probably as much as 10 times higher than those 
e no experienced in the real world. As the discussion of hi-valu results will indicate, this bias had 
de- to be considered and used in the interpretation of results. This is because the hi-valu defer- 
ant was ment policy was designed to take account of the fact that most hi-valu parts have very low 
demand rates. The occurrence of unexpectedly high demand rates caused all sorts of diffi- 
ndi- culties in the operation of the hi-valu policies for LS-2. 
any of In the Category 2 and 3 parts, a sample was used which had too few zero-demand 
em, items in comparison with real-world experience. Although this is of some significance, it is 
under - not believed to have had important effects on the Category 2 and 3 results. The main effect 
data is that we could not evaluate the results of stocking the zero-demand items as the LS-2 Cate- 


gory 2 and 3 policies would tend to encourage, especially if the part was cheap enough. 
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Because biased demand rates were used for both the hi-valu and Category 2 and 3 





parts, it is difficult to determine if the real-world experience of the floor personnel was used 
as effectively as it might have been had the data used in the simulation been more realistic. 
We do know that the floor managers stated many times that the consumption rates they were 
experiencing in LP-I were different than those they would have expected in the real world for 


the same kinds of line items. 


Factory Simulation 





The factory simulation used in LP-I was operated manually because not enough time 
was available to program the calculations of the factory for the electronic computer. This 
meant that many errors in record keeping and computing necessarily crept in because of the 
thousands of manual calculations and entries that had to be made with the factory records. 

There also were policy failings in the factory simulation, especially for hi-valu parts. 
For example, the factory did not afford the same degree of priority to ANFE!” items as it did 
to AOCP!” items. Although the designers of the hi-valu policy in LS-2 had thought that ANFE 
items would be treated with the same importance as AOCP items, they had not made this an | 
explicit part of their policy statement, and so the factory felt it was not bound to that interpre- 
tation. Also, the factory misunderstood the super-expedite!® arrangement which had been made { 
with the LS-2 hi-valu manager, and it tended to deliver hi-valu parts on a super-expediate basis 1 


much more slowly than the management personnel had thought it would. 


Personnel Biases 





Undoubtedly many personnel biases also crept into the LP-I simulation. For example, 
there was no replication or duplication of the personnel manning the several important 
management positions. It is therefore difficult to determine how much of the effect of a policy 
was due to the policy itself and how much was due to the individual who was using that policy. 
If more opportunity had been available to test the policies with different individuals, thensome 
distinction between the effect of humans and the policies in LP-I could have been made. This 
is especially true of the provisioning conferences because many of the decisions made in 
these conferences determined much of the outcome of the subsequent experimental experience. 

Other human errors also affected the experimental results, including errors in 
interpretation of procedures, errors in use of data, and errors in simulation. Also, because 
so many human beings were used in the critical positions in the laboratory, and they were so 
active under compressed time, it was not possible for the laboratory staff to be completely 
sure that the personnel in the critical positions conformed to policy. Certainly they had the 
opportunity to avoid or to violate policy if their feelings went in that direction. However, we 
cannot be sure that their violations in LP-I were any more frequent or extreme than those 


they commit in the real world. 
I7ANFE = aircraft not fully. equipped; lacks parts that affect combat performance of aircraft, 
such as the fire control system, but do not prevent flight. 
AOCP = aircraft out of commission for parts (parts preventing flight). 19 
8super- expedite means that the factory bent every effort, short of disrupting the aircraft 
production line, to supply spare parts. 
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Single Run 


sed Also, LP-I contained a single complete run. Thus, the comparison between LS-2 and 

a LS-1 in over-all terms, such as in cost and effectiveness, is based on a sample of one. If 

re different people and a different sample of input data had been used, the results might have been 

for somewhat different, but how different is a moot point. At the same time it should be recog- 
nized that, in terms of specific decisions resulting from use of the supply policies, the mana- 
gers produced large samples of experience. Thus, the managers applied the hi-valu policies 

- toa sample of approximately 50 line items over a 3.5 year period, and they applied the 

. Category 2 and 3 policies to over 700 line itemsfor a similar period. In this sense the sample 

_ size of separate procurement and distribution decisions amounted to several thousand elements. 
Under these circumstances, the over-all results are probably quite reliable, even though the 

ste. experiment contained a single run of both logistics systems. 

it did Results of Category 2 and 3 Policies 

NFE Table 1 shows the expenditures made by LS-1 and LS-2 for Category 2 and Category 3 

an parts. It is seen that LS-2 spent slightly less than LS-1 for such parts. The table also shows 

rpre- the number of base stockouts!? suffered by LS-1 and LS-2 for Category 2 and 3 parts. Here, 

n made too, in the case of these actions, LS-2 showed to far better advantage than LS-1. Therefore, 


te basis we can say that LS-2 did significantly better than LS-1 in the support of Category 2 and 3 
parts. This was shown by LS-2 spending less money and achieving far greater support effec- 


tiveness than LS-1. 


imple, TABLE 1 
Category 2 and 3 results 
(14 Quarters) 





























policy 

licy. Expenditure 

anes A+tion (in thousands of dollars) 

This LS-1 LS-2 

l Category 2 

“ience. Base stockouts 4,909 2,955 
Spares expenditures 1, 668 1,491 

cause Priority requisitions 6,562 746 

re so Routine resupply actions 7,089 1, 960 

tely re ——_---————————_}——_ 

i the Category 3 

r, we Base stockouts 3, 121 544 

ose Spares expenditures 262 240 
Priority requisitions 7,100 285 

craft a 2. = 

- '9A stock-out occurred when a part called for in the repair of an aircraft could not be provided 


‘rom the supply or maintenance on the base. 



























M. A. GEISLER 





Table 2 illustrates the effects of the different policies used by LS-1 and LS-2 admin- 
istering Category 2 and 3 parts. As can be seen, the reorder points and stock control levels 
of Category 2 and 3 parts were much higher in LS-2 than in LS-1. These higher levels in 
LS-2 had the effect of producing much more protection for the bases in meeting demands for 
Category 2 and 3 parts and in reducing the number of times the bases had to be resupplied with 


Category 2 and 3 parts, either routinely or under priority conditions. 


TABLE 2 
Base Stockage Policies in Category 2 and 3 























Average Quantity Over Selected Line Items 
Stockage Catetory 2 Items Category 3 Items 
Factor 

LS-1 LS-2 LS-1 LS-2 

Reorder Point 1.0 5.1 0.6 4.6 

Reorder Quantity 0.6 5.6 1.8 11.8 
Stock Control | 
Level 1.6 10.7 2.4 16.4 | | 

















Results of Hi-Valu Policies 

The upward bias in the hi-valu parts previously described had a considerable effect on 
the results obtained in the LS-2 use of the hi-valu policies. Table 3 shows the cumulative 
percentage of items in the real world (Stock Balance and Consumption Report) and the corres- 
ponding cumulative percentages in the LP-I sample by demand-class interval.29 We have 
also shown the expenditures and the stockout days of LS-1 and LS-2 for the line items con- 
tained in each of these demand-class intervals. By examining the behavior of LS-1 and LS-2 

TABLE 3 


Hi-Valu Results — Spares Expenditures And Stockouts 
(14 Quarters) 




















Number of : ; 
Demands | immer | CUGr tine ome: | (Thousands of Dallas) | Scout Days | |) | 
| Flying Hours | Sample) ppy | sB&CR|  LS-1 ts-2 | Ls-1 | Ls-2 L 
0-1 - 15.1 60.1 | 71 21 4 7 | 
1.1-12.5 9 32.1 90.2 | 432 72 | 95 108 
12.6-25.0 © 47.2 94.8 | 260 361 239 243 
25.1-50.0 5 56.6 97.2 | 119 151 68 102 : 
50.1-100.0 10 75.5 99.4 | 292 426 553 331 
100 and over 13 100.0 100.0 | 606 943 886 | 1369 _ 


























20Since real parts were used in LP-I, it was feasible to go to the Stock Balance and Consump- 
tion Report and establish the percentage distribution of the sample in the different demand- th 
class intervals. It was thought that these were the data being used until the bias was detected 
and the reason for the bias learned. 
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by demand-class interval, we can try to estimate how LS-2 would have functioned if the de- 

































































“ mand rates had not been biased. 

As can be seen, LS-2 spent less than LS-1 on hi-valu parts up to a demand rate of 25 
per 100,000 flying hours. Beyond this demand rate, LS-2 spent more than LS-1; the extra 
wet expenditures of LS-2 were largely for priority actions, such as expedited production and trans- 

portation, caused by larger demand rates than those expected. 

According to the Air Force Stock Balance and Consumption Report, 95 percent of the 
hi-valu line items in the real world would have had demand rates of less than 25 per 100,000 
flying hours. Consequently, if the LP-I sample had been unbiased, LS-2 probably would have 
spent less on 95 percent of the line items. 

At the same time, LS-2 had about as many stockouts as those of LS-1 for the line items 
with demand rates under 25 per 100,000 flying hours. Therefore, we can say that the hi-valu 
policies worked well for parts that had demand rates of 25 per 100,000 flying hours or less and 
that in the real world this applies to the bulk of hi-valu items. This result accords with 
expectation because it was realized that beyond a particular demand rate deferment of procure- 
ment could not be used, and for these high-demand parts a system more like that used current- 
ly in the Air Force would be appropriate. Because comparatively few parts have such high 
demand rates, valid estimates of their demand rates can be obtained quickly, and, therefore, 
for these few parts, surpluses and bad procurement should not be a problem.2! 

= Table 4 shows the initial provisioning experience of LS-1 and LS-2 on the C and D 
, series of the aircraft used in LP-I, and it highlights the differences in the hi-valu policies 
ial used by LS-1 and LS-2. As the table indicates, LS-2 spent very little money in its initial 
ad provisioning conferences, whereas LS-1 spent large sums of money, sums which were the 
a , bulk of expenditures for spares of hi-valu items. 
TABLE 4 
Obligations Made During Provisioning Conference for Hi-Valu Items 
Number of 
—_——, | Demands Cumulative Percent "C" Series of A/C "D" Series of A/C 
Days ad aa 000 of eine. ems LE (Thousands of Dellers} (Thovaente of Deters 
| _ Flying Hours I LS-1 LS-2 LS-1 LS-2 
LS-2 i 0-1 60.1 | 15.1 19 5 14 1 
a | 1.1-12.5 90.2 32.1 213 sj 84 + 
1 | _ 12.6-25.0 94.8 | 47.2 29 9 81 4 
108 |__25.1-50.0 98.2 | 56.6 52 3 54 2 
243 50.1-100.0 98.9 | 75.5 65 4 78 7 
102 | Over 100.0 100.0 | 100.0 174 34 224 0 
331 Sa Piraae 
1369 
ial al, simple weighting scheme was used to estimate what the aggregate expenditures for LS-1 


and LS-2 would have been had unbiased demand rates been used. This weighting indicated 

and - that LS-1 would have spent about 900 thousand dollars and LS-2 about 500 thousand dollars, 
tected which represents a considerable difference in favor of LS-2; in addition, the number of stock- 
outs would have been about the same in both logistics systems. 
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Data Processing Center Results 

LP-I experience showed that the use of a data-processing center is a basically sound 
concept. The DPC fitted in well with the new supply policies for hi-valu and Category 2 and 3 
parts used by LS-2. No significant amount of conflict was noted between the base managers 
and the weapon system managers in LS-2, despite the fact that the bases had no record of 
current or due-in assets but depended on automatic resupply from the storage site. However, 
it did take time for the base managers to become accustomed to this new situation, and there 
was considerable discussion during the early phases of LP-I as to whether or not the bases 
could operate realistically without such records. 

Also, the complete LS-2 system, including the DPC operation and the centralized com- 
putation of levels and requirements, was very difficult for the floor personnel to understand. 
As a result, much Statistical Service activity was required, both to explain how the system 
operated and to investigate apparent discrepancies which were reported by the floor personnel. 
Many of these discrepancies were not errors but were unexpected outcomes which the floor 
personnel had difficulty in reconciling with their previous experience. 

It is clear that the use of a system like LS-2, with its new policies and the data- 
processing center, will bring with it much Statistical Service activity until the many personnel 
who have contact with the new system feel at home with it. Its centralization of computations 
makes the new system a more delicate one, especially since some of the policies depend on 
the responsive character of electronic computations. This means that there will have to be 
continual quality control in the operation of the system, and trouble spotters will have to be 
employed to insure that the system continues to function properly if the full benefits of this 


new technology are to be realized. 


Implementation Experience 


Category 2 and 3 Policy 





Much implementation experience in the use of the new policies was obtained in LS-2. 
Several examples of such experience will now be cited, not so much to emphasize their sub- 
stantive meaning but to indicate the important procedural detail that was learned in the LP-I 
type of logistics simulation. For example, in the use of Category 2 and 3 policies, it became 
clear that more allowance for differences in logistics requirements among bases must be made 
in the formulas computed at the data-processing center. These differences should take account 
of the variations in demand rates, geography, and field-maintenance practices. 

Also, changes should be made in the way demand rates for the parts repair depot and 
the IRAN activity were obtained in LP-I. These activites should be treated as separate pro- 
gram elements so that separate consumption rates and program data for each activity are used. 
In LP-I a single over-all consumption factor and program element were used. 

It was also noted that there was a significant difference between LS-1 and LS-2 in the 
amount of amended shipping instructions““ for Category 2 and 3 items. LS-1 had only 25 


22an amended shipping instruction directs the factory to change the previously established 
destination for a shipment. 
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amended instructions, while LS-2 had about 500; thus a large administrative workload was 


und imposed on the LS-2 factory manager, a fact which should be recognized. This came about 

ind 3 | because LS-2 used the factory as a sort of alternate storage site. When an order was first 

rs placed, the delivery address given was thé storage site. When the order was ready, the de- 

f cision was made as to the best destination for it. Often, this resulted in a change in the 

ever, delivery address. Usually, the order was delivered to a base, since the LS-2 Category 2 and 

ere 3 policies kept very little inventory at the storage site but rather placed almost all of the in- 

eS ventory in the bases. 

aise Hi-Valu Policies 

—y Important implementation experience was also obtained with hi-valu policies. For 

-_ example, the LP-I experience strongly indicated that a differential buffer policy should be 

nui used, which takes account of the differences in characteristics of the individual hi-valu items. 

_ In LP-I, a single-standard buffer policy was applied indiscriminately to all the hi-valu items. 
A better understanding of the differences between insurance and operating items was 

aes also obtained because of the LP-I experience. ” Also, better rules for deciding whether an 

— item is an insurance item or an operating item were obtained, and better techniques were de- 

Mean veloped for reviewing the experience of each line item to determine its type. 

bitte The sensitivity of the hi-valu policies to due-in information was also understood better 

set because of LP-I experience, and the importance of receiving accurate and timely due-in in- 

stoi formation from the factory system was found to be a most important requirement of these hi- 

in valu policies. 

Also, the LP-I experience indicated that the factory must have a flexible ability to fill 
orders for hi-valu items with little or no advance notice. The uncertainty as to when a parti- 
cular order would be placed by the supply manager makes it necessary for the factory to have 
this flexibility. Also, the manager must have the option to change procurement schedules fre- 

S-2. quently as experience unfolds. 

sub- The experiment also indicated that improvements are desirable in the way IRAN re- 
LP-I quirements for hi-valu items were obtained in LP-I. It was indicated that rather than charg- 
ecame ing the consumption of parts used in IRAN back to the base from which the aircraft came, it 


ye made would be better to maintain a separate IRAN consumption record and to use the number of air- 


account craft entering IRAN as a program element in the computation of hi-valu item requirements. 
LP-I 

rt and wan 

, pro- | One of the important functions of the laboratory is to test proposed policies under spe- 

a ° cial conditions and circumstances that cannot readily be duplicated in the real world. The 


testing of new policies under wartime conditions was one of these special contributions intended 
from the laboratory during LP-I. 


in the 

25 23an insurance-type item is one with an expected or experienced very low demand rate. An 
operating item is one with an expected or experienced relatively high demand rate. The hi- 

hed valu policies operated differently depending on whether the item was an insurance or oper- 


ating item. Important procurement decisions hinged on this distinction. 
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Design of LP-I Wars 

During the LP-I exercise the logistics systems managers were directed to plan for war 
as they would do generally in the real world. They were given simulated war planning data and 
other guidance to use in establishing a war readiness comparable to that required by the real 
world policies. Both logistics systems were given the same war planning guidance for their 
use during the LP-I run. 

To facilitate the laboratory operation and to obtain valid results, both for the peacetime 
activity of LP-I and the wartime activity, the design selected was to have both logistics 
systems run through the peacetime program that was to be simulated in the laboratory. In the 
meantime, two dates were selected as D-days for each of the two wars which would be later 
simulated. After the peacetime run was over, two separate wars were run. The first war 
began without warning on the first day of Year 2. The second war began, with 5 days of warn- 
ing, at the start of Year 3. We felt that two possible ways in which war might begin should be 
tried to get a better picture of how the policies would work in wartime. 

The simulated wars were assumed to be thermonuclear. Under these circumstances, 
certain of the bases simulated in LP-I were assumed to be destroyed through attack. In addi- 
tion, it was assumed that support from the depot was not teas.ole during this period, which 
required putting stress on local base resources and improvisation. At the same time, the 
transportation and communications systems were assumed disrupted so that normal use of 
these facilities was denied to the logistics systems. It was assumed that the data processing 
center used in LS-2 continued to operate. 

Each of the two wars lasted for an arbitrary 10 days. On each of these 10 days the 
surviving bases flew their aircraft at maximum sortie rates. Each base could send its in- 
commission aircraft on sorties four times daily: at 0800, 1000, 1300, and 1500 hours. 


24 was assumed 


The logistics situation was as follows: first, organizational maintenance 
to continue around the clock; second,*> field maintenance was increased by 50 percent over its 
peacetime level, mostly by the use of overtime; third, the failure rates were changed to reflect 
wartime utilization of aircraft; fourth, the bases were assumed to be dependent on their own 
stock, on cannibalization, or on borrowing from other bases, since depot support was not 


available. 


LP-I War Results 
The general results of both wars were that both logistics systems performed equally 

well. Neither system seemed to fly significantly more sorties than the other system, although 
LS-2 did seem to be able to meet its logistics requirements from the D-day serviceable in- 
ventory, with less recourse to cannibalization, use of next higher assembly, exchange of parts 
between bases, or intensive use of field maintenance. 

240rganizational maintenance is that maintenance performed directly on aircraft to make 
them combat-ready. 


5Field maintenance is that maintenance required to restore unserviceable spare parts to 
serviceability. 
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It may be of interest to go a bit deeper into the specific results of War 1, since it 


































ee: represented a more severe test for the hi-valu policies. The deferment practiced by the hi- 

a and valu policies of LS-2 was most marked in the early periods of the phase-in, so that it was of 

eal interest to see how LS-2 functioned at the end of the first year, when War 1 occurred. It was 

sir found that both systems flew roughly the same number of sorties, although LS-2 did slightly 
better. Also, LS-1 suffered many more stockouts than LS-2. Breaking down the stockouts by 

etime price of the line items affected, we find that LS-2 had more stockouts with hi-valu parts, but 
LS-1 had so many more stockouts with Category 2 and 3 parts than LS-2 that the over-all 

n the number of stockouts of LS-2 were significantly less. 

ad The LS-2 bases were especially better stocked with Category 2 and 3 parts, for which 

ad the bulk of the demand occurred. LS-1 was much more dependent on cannibalization and field 

_— repair of spare parts because of its more limited base serviceable inventory for meeting the 

ld be parts demands that arose during the war. 

sini Concluding Remarks 

addi- The results obtained in the simulation of both wars seemed to be quite reasonable. 

ch This was the first significant attempt by the laboratory to simulate war, and it can therefore 
be considered as a maiden attempt. Undoubtedly, future war-simulation efforts can be much 

. more sophisticated. In LP-I for example, we might have varied the daily sortie activity re- 

; quired of each logistics system. It is not necessarily true that interceptor aircraft bases 

— would be required to fly maximum sorties each of the days of the wartime period, as was 
assumed in LP-I. No allowance was made in these wars for partial damage, and this is also 

: likely to be characteristic of real war. Finally, the assumption of no depot support during 


the war period may not have been completely realistic. 


ssumed LABORATORY MANAGEMENT EXPERIENCE 


ver its Close Tie Between Policy Research and Simulation Activity 

reflect The following brief remarks are intended to provide some indication of what was 
own learned during the course of LP-I in the management of a laboratory such as the Logistics 
yt Systems Laboratory. These remarks might give the reader some feeling for the complexi- 


ties involved. A most important experience was to recognize the great need for close com- 
munication between the research activity which specifies the policies to be tested or studied 


in the laboratory and the laboratory simulation activity. There must be a very close com- 


ally munication between these two efforts to insure that the simulation faithfully studies what 
ithough should be examined with the policies and that data which are gathered during the test will be 
in- of maximum usefulness to the policy makers in revealing to them the ways in which their 

f parts policies may work better in the real world. 


Required Flexibility of Laboratory Organization 
i The laboratory organization must be very flexible because it must change continually 
during an experiment. In the early stages of the experiment's development, the organization 


is oriented very much towards the development of the mathematical models that form 
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the basis for the simulation and the gathering and development of the data necessary as inputs 
for the experiment. This period also involves considerable design of analysis systems and of 
other means by which the laboratory can extract requisite information, both during and after 
the experimental runs. After the simulation is designed, the organization turns to the more 
operational procedures, including the development of the manuals that guide the participants, 
the preparation of training courses, and the rehearsing of the simulated system for flaws be- 
fore the actual laboratory run begins. Finally, when the run actually starts, the laboratory 
organization is concerned with operation of the experiment, with collection of data, and with 
analysis of data. 


Scheduling of Experiments 

The scheduling of a laboratory experiment must be very realistic, and it must take into 
account the amount of time required to prepare each stage. In particular, this includes 
allowing enough time for the programming of the electronic computers and the rehearsal of the 
system before the floor personnel are permitted to operate the system. A study of a complex 
logistics system involves so many interrelated and important elements that, unless there is 
adequate testing and rehearsal of the system, many flaws might go undetected, and they could 
persist throughout the experiment, with resulting loss in benefits and results obtained from 


the experiment. 


Short Runs 

The LP-I run lasted more than 3 months. The experience during LP-I indicates that 
the runs for each experiment should be much shorter, so as to sustain greater interest on the 
part of the participants in the experiment and to reduce the laboratory's reliance on a single 


long run that might contain serious errors which could ruin an experiment. On the whole, 


several shorter runs should produce a larger net payoff for laboratory experiments. : 

Desirability of Experiments 
The type of laboratory experimentation used in LP-I is a necessarily expensive under- - 

taking. However, it is the only presently known methodology for studying the kinds of problems 

that are so important in a logistics system. These problems are the implementation of com- 

plex policies that affect many organizations in the logistics system and the study of the organ- . 

izational interactions, data flows, and problems that are raised through the use of these S 

policies. The lessons learned in a laboratory will do much to avoid wasteful operation in the 

real world, where the costs and consequences are much greater than in the laboratory. 

Analysis ' 
A tremendous amount of data is generated during an experiment such as LP-I. This ( 


means that ample provision must be made for an adequate analysis of these data. It is well 
and desirable that the personnel from the simulated organization (e.g., Air Force personnel 
in LP-I) participate in such analyses so that they can better realize the intentions of the 
experiment and how well these intentions were accomplished during the run. Their familiarity 
with the experiment permits such personnel to interpret and understand the results for 
application to their real-world organizations, which is an important part of the necessary 


bridge between research and implementation. 





* * * 

























AN OPTIMAL INVENTORY SOLUTION FOR 





nts, SOME SPECIFIC DEMAND DISTRIBUTIONS* 
be- 
ry Brian Gluss 
ith Armour Research Foundation of Illinois Institute of Technology 
nto = ee ee OO | 


A multistage inventory decision process is considered in which the 
penalty cost for failing to have in stock sufficient items to satisfy de- 
of the | mand at the next stage includes a fixed administrative cost and a cost 
nplex proportional to the deficiency in stock. When the demand probability 
density function belongs to a certain class, which includes all Pearson 


dos distributions, a method is found for determining whether or not a con- 
could stant stock level solution is optimal. 
om L : = es 
hat 
or INTRODUCTION 
angie The problem has been considered [1] of optimizing an infinite-stage decision process 
in which stock is ordered at the beginning of successive periods of time to meet demands 
whose probability density functions 9(s) ds are known. If the demand proves to be greater 
than the stock x in hand, a "penalty cost" is incurred in rushing through an order to meet the 
excess demand, whereas if the demand is very small relative to the stock in hand the process 
der - might again be uneconomical. Bellman has shown that, in the case where the cost of ordering 
omen z items to increase the stock level is kz, and the penalty cost of rushing through z items to 
ee meet a demand of z over and above the stock level is of the form pz +q, the optimal 
. solution consists of the two regions, (0, X) and (xX, ©), in which the ordering policies are: 
. the (i) Increase the stock to X O<x<x 
(ii) Order nothing X > X, 
provided that there exists only one minimum, at the point x, of the function; 
{ ¢-@ eo y ) 
me (1) Fly)=ky+a ip f (s-y) 9(s) ds + | £(0) + q] j 9(s) ds uf f(y-s) 9(s) ds}, 
vell L “y y 0 J 
ynnel 
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where a, (O<a< 1), is a fixed discount ratio used in each period for discounting future costs 
and f(x) is the expected total discounted cost when one starts with an initial supply x and uses 
an optimal policy. 

In [ 1] a sufficient condition, ¢(s) < 0 for all s, is given for a solution of this kind. In 
the present paper, less stringent requirements are made upon 9¢(s), and a simple procedure 
is obtained for determining whether or not Pearson-type (see [ 2]) demand functions 9(s)ds 
give the above policy solution. 


THE PROCEDURE 


All the relative minima of F(y) satisfy the equation 


(2) G(y) = (1 - a)k - a(p - k) 7 o(s)ds - aqg(y) = 0. (See [1] .) 
The roots of this equation are vena by those of 

(3) Gy = MP - k)Gy) - agg"(y) = 0; 

that is, 

(4) aa — xk. A, say. 


Let it be assumed that G(0) < 0, that is, ap >k - aq@(0). This is not a very restrictive 
condition and is satisfied in particular by all cases where the discounted penalty cost per item 
ap is greater than the initial ordering cost per item k. 

From (4) it is seen that whether or not there exists only one minimum of F(y) is 
easily determined for any demand function 9(y) of the form 


sy) LG) 

dy) MG)’ 

where L(y), M(y) are quadratic functions of y. When L(y) is linear, 9(y) is Pearsonian, by 
definition, and the results of this paper therefore apply to such distributions as the Normal, 


t (Pearson type VII), x2 (Pearson type III), and F (Pearson types depending upon the number 
of degrees of freedom). 


(5) 


The procedure that is followed is to find the real positive roots, if any, of 
(6) L(y) = A M(y). 


G(0) = - (ap - k) - aqg(0) <0, and G(~) = (1 - a)k > 0. Hence if there is at most one positive 
root Yj of (6), there is only one solution to (2) and, depending upon whether G'(0) = 


a(p - k)¢(0) - aqg'(0) is greater than zero or less than zero, the solution to (2) lies in the inter- 
val (0, y,) or (yj, ©), respectively. (See Figure 1.) 


Figure 1 
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If there are two positive roots y;, yg (yy <yg) of (6), however, the number of solutions 
to (2) depends upon the sign of G'(0): 

(i) G'(0) <0. 

In this case there is one solution of Eq. (2), and this lies between y; and yo, as can 


be seen from Figure 2. 


| ae 
| - (1- adi 
Figure 2 + >~y 


(ii) G'(0) > 0. 


Here the number of solutions of Eq. (2) is one or three according as G(yo) is respec- 





tively greater than or less than zero. (See Figure 3.) 








To summarize, there is one root of Eq. (2), and the results of this paper therefore 
apply, when 

(a) there is one positive root of Eq. (6), 

(b) there are two positive roots of Eq. (6) and G'(0) <0, or 

(c) there are two positive roots of (6), G'(0) > 0 and G(yg) > 0. 


Two examples are now given to demonstrate the procedure. 


Example 1: The Chi-Square Distribution 
By definition, 
oy) —— “s @ty<. 
2 
It will be assumed that n > 2. Hence, 


oy) _n-2-y 
oy) 2y : 





and the solution to (5) gives one root (n - 2)/(2A + 1). If follows that there is a con- 
stant stock level solution, with x in (0, y,); x may then easily be found from (2) by successive 
approximation. 

For example, consider n= 12, p = 2k, q= 10k, anda=0.97. Then A= 1/10, and 


y; = 0.8333. By successive approximation the value x = 7.848 is obtained. 
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Example 2: The F-Ratio Distribution 








" 
¥ . 7 
0< < 
Hy) @ ™ n+ Rp <y 
1 + To y 2 
n n 
M1) ng -( (744) my 
Hence, oy) \2 2 
9(y) y(ng + nyy) I 
and (6) reduces to of 
n n _ 
(7) Anyy" - Ana( +1 Jo y ~ ) = ©. 


It follows that, for ny >2, there is at most one positive root y,, and hence the of 
solution is a constant stock level solution. et 
ti 
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DETERIORATION OF INVENTORY AND EQUIPMENT* 


M. Klein and L. Rosenberg 


Columbia University 


INTRODUCTION 
The purpose of this survey paper is to describe problems relative to the management 


! It is traditional [16] to dichotomize this class of items into a 


of goods which deteriorate. 
"constant efficiency"' and a ''diminishing efficiency" group. 

Constant efficiency goods are characterized by the stability of the unit value of their 
output over their lifetimes. The deterioration is basically expressed in terms of the lengths 
of their lives; these are usually taken to be random variables. Typical constant efficiency 
goods are electric light bulbs and fuses. 

A diminishing efficiency item is one whose net output value decreases with use or with 
time. The deterioration of output or efficiency is caused by wear, maintenance costs, obso- 
lescence, etc. The life of goods in this class can usually be extended almost indefinitely by 
appropriate maintenance measures. Most industrial equipment as well as many kinds of 
inventory items fall into this category. Actually, there are goods whose value may increase 
with time in storage, e.g. whiskies, works of art, etc. Problems involving goods of this 
type also have been treated. 

The problems dealing with equipment, discussed below, are scheduling problems. 
They are concerned with the optimal timing of decisions to replace or repair a piece of equip- 
ment which deteriorates in some known way (problems 1 and 2) or in some probabilistic way 
(problems 3-5). 

Work done on the management of deteriorating inventory also is discussed. The prob- 
lems are of two kinds: one involves the determination of optimal policies with which to de- 
plete inventories, the other is directed toward finding sampling inspection policies to aid in 
the maintenance of prescribed quality levels of stored goods. 

There is a large body of literature which deals with the problem of estimating the de- 
terioration function of either constant or diminishing efficiency goods. This paper will not 


review this material; however, an excellent bibliography | 17] is available. 





*Manuscript received May 7, 1959. Research was supported by the Office of Naval Research 
under task NR-042-034. 

Most of the problems included have been considered in the literature; problems 3 and 4 do 
not seem to have been stated previously in the forms given here. 
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Reliability problems usually involve some form of deterioration; hence, they could 





properly have been included in this survey. However, to keep the task manageable, we have 
included only some problems of the reliability type which involve single components. A larger 


mn 


class, concerned with the management of multicomponent systems, is not discussed. 


EQUIPMENT DETERIORATION 
This section describes some problems which arise when equipment deteriorates in use. ( 
Our prime interest is in finding optimal replacement schedules for such equipment. Some- 
times it is possible to repair a piece of equipment or renovate a stored item so that it can be . 
returned to an "'as new" condition with respect to all cost and performance parameters. In P 
this sense, the models used in problems 3-5, described below, may be considered as suit- 
able for both equipment and inventory management. On the other hand, the models discussed 
in problems 1 and 2, in which there are decreasing salvage values and increasing maintenance 
costs, reflect what may be termed irreversible deterioration properties; these commonly 


apply only to equipment problems. Some inventory management problems involving storaged 





goods which deteriorate in this irreversible fashion are discussed later in this paper. b 

Problem 12 P 

As an introduction, we start with a description of an elementary form of a classical 
equipment replacement problem. A piece of equipment whose performance is known to deter- P 
iorate with use is to be put into service. The deterioration is assumed to be entirely reflected " 
in the operating costs. Suppose that the cost of the equipment is denoted by A and that the rate . 
at which operating costs are incurred can be represented by a continuous function of time m(t). . 
It is of interest to determine the time (t,)) at which a new piece of equipment with identical a 
characteristics should be put into service. This value ty is roughly what is usually called the “ 
"economic life" of a piece of equipment. 

The average cost per unit time associated with a piece which is replaced at the end of 
its economic life, to, is given by 

Ao) = +E ‘y m(t) dt. 
O 
If m(t) is nondecreasing, the (unique) minimum of Q(to) occurs when 
m(tg) = = + c (° m(t) dt. 
fe) 

Hence, the optimal replacement time occurs when the operating cost rate is equal to the 
average cost. 5 

Sometimes situations occur in which the operating costs are linear [5]. For this case, a; 
if m(t) = bt, the economic life takes the simple form ra 

VD e 





2See [ 6] for a discrete time version of the problem. 
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More generally, an equipment replacement model includes consideration of two other 


ave components in which deterioration plays an important part. These are the salvage values, 
arger s(t), and the revenue from the equipment, v(t). Thus, if it is assumed that no technological 
changes will occur, the problem is to find that value of tp which maximizes the present value 3 
of the total profit, 
n. == t 
oma (1)  G{to) = a e Kto a + J [v(t) - m(t)] —— dt - S(t) _ tte ( 
1e - 
> tes where r denotes the rate of interest and n denotes the number of replacement periods in the 
In planning horizon. 4 The presence of the discounting factor e “rt permits consideration of an 
it- infinite planning horizon, i.e., n= o. 
—- An interesting characteristic of this problem is the reported insensitivity> of the cost 
on function (1) to changes in the replacement interval around the region of the optimal replacement 
y time. 
aged Another approach to the determination of an optimal equipment replacement policy has 
been suggested by Bellman (2) and Dreyfus (10). Their method is especially valuable when 
predictable technological changes are anticipated. 
_ Consider a requirement for a policy for a finite planning horizon consisting of N equal 
Rates periods. At the start of any period the decision maker is faced with two possible choices: he 
flected may decide to replace the equipment with a new piece (decision R) or keep the old piece (deci- 
—* sion K). Clearly, an optimal replacement policy (a series of N decisions) must have the prop- 
- erty that at any time (n), no matter what policy has been followed in the past, the remaining 
— decisions must be optimal. That is, at any time it is desirable to choose the alternative which 
ed the maximizes the return for the remaining periods in the planning horizon. Thus, let 
t = the age of the equipment in use at the start of the period under consideration 
ee fy(t) = the profit, associated with equipment of age t at the start of the nth period, which 
is obtained during the remaining periods in the planning horizon (n,..., N) if an 
optimal policy is used for these periods. 
vy(t) = revenue during the nth period from a machine of age t. 
m,(t) = maintenance costs during the nth period for equipment of age t. 
c,(t) = cost of replacing the equipment at the start of the nth period if the old equipment is 
of age t(i.e., the purchase cost minus the salvage value. 
1e 
—_— "ht preciect value of a profit made at some future time is the amount which one has to set 


aside at the present, drawing compound interest at the rate r, which will be equal to that 
future profit. (See [ 3] p. 315.) 

For further discussion of this model and others in the same spirit see Hotelling [13], Roos 
(20) , and Preinreich [19] . See Smith [22] for a numerical example. 

See Clapham [5]. Smith [22] and Bowman and Fetter [3] . 
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Then, an optimal policy at the start of the nth period with equipment of age t has a 
present value equal to 


R:vp(0) - mp(0) - ep(t) + af, | WF, 
fy(t) = max 


R,K K:vp(t) - mp(t) + a fp 41(t+1) 


where t= 1, 2,..., n anda denotes the appropriate discount factor, An optimal policy may now | 


be obtained by using this recurrence relationship successively through the N periods, starting 


at the last period and noting that a finite planning horizon implies that fy, ;(t) = 0. The sequence 2 


of choices giving rise to f, (t*), where t* denotes the age of the equipment at the start of the plan-)7 


ning horizon, is the optimal policy. In general, it will not involve periodic replacement. How- 
ever, for the case of no technological change, the functions v;, my, and cy do not change with 
time, and the optimal replacement policy takes the form of periodic replacement. 

In addition to the conceptual and computational simplicities offered by the dynamic 
programming formulation of the problem, there is an advantageous flexibility inherent in this 
kind of multistage structure. At each stage, if desired, additional alternative decisions may 
easily be introduced. Thus, one can incorporate a series of different maintenance policies 
as subalternatives under decision K. 

A troublesome problem associated with the use of any of these models is the require- 
ment for somewhat extensive information concerning the future behavior (deterioration) of the 
equipment and of the market in terms of the functions, v, m, and S. As an example of diffi- 
culties involved, the determination of v requires knowledge of the output of the equipment, of 
the quality of the output, and of the value or price of the product over time; the prediction 
of the output in turn usually calls for a prediction of demand. If the equipment has been used 
in the past, estimates of the maintenance costs and salvage value may not be a serious prob- 
lem. However, the prediction of these latter costs may be quite irksome where there has 
been no previous experience with the type or model of equipment of interest. 

Usually, the use of an optimal maintenance policy is assumed when models such as 
(1) are employed. One approach to the determination of an optimal intraperiod maintenance 
policy is discussed next. It may be noted that the model used below is actually a special 
case of (1). 

Problem 2° 

An item deteriorates in use or in storage until it reaches a point at which it can no longe! 
perform its intended function. Before it reaches this point it can instantaneously be serviced to 


return its condition to that equivalent to a new item. The cost of such service is assumed to be 


a function of its condition (or, equivalently, of the time since its last servicing took place). Afte 


a certain period it becomes obsolete and is replaced, regardless of its condition. 


6This problem has the same structure and solution as one of Savage's ''cycling problems'' 
[21]. 
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The problem is to determine the times at which the item is to be repaired subject to the 








a 
restriction that, except for the servicing time, the item is to be kept in operating condition. 
Let xj denote the time between the (i-1)th and the ith servicing, take the total interval 
of interest to be unity, and let n be the total number of servicing times, i.e., X; is the 
i= wi 
time at which the last servicing takes place. Since the item is kept in operating condition, 
there exists an upper bound, denoted by b, to the length of time between servicing intervals. 
We will also assume that f(x), the cost of servicing after time interval x, is a nondecreasing 
may now 
function. 
sate. The problem is to determine the values} x, t j(i=1, 2,..., n), and "n" which minimizes 
rs the total cost, 
the plan- 
How- (2) Q(x}; og aes n) > > f(xj), 
i=1 
‘e with 
subject to the following restrictions: 
ic yx =1 
n this = 
} may 0<xj <b 
-ies If f(x) is differentiable and known, the problem can be solved by the usual methods of 
calculus to enumerate optimal solutions to (2) for each of the possible values of "'n". 
juire- Some slightly better results are obtainable for the following special conditions: 
of the Case I: f(x)is convex. Here it is easy to show that the minimum is obtained for fixed 
diffi- "n" by taking all of the servicing intervals to be equal; i.e. xj= 1/n, i= 1, 2,..., nfor 1/n<b. 
nt, of Hence, (2) reduces to 
1 
a (2") Q(xy, +--+, Xp) =nf 2 
1 used 
prob- and n may be computed as one of the integers which minimizes (2'). 
_— Case II: f(x) is concave. Since Q will also be concave, the minimizing solution occurs 
at an extreme point of its domain. Hence the optimal program is to have the intervals between 
_— servicing as long as possible. Thus, we choose x; = b for as many intervals as necessary, and 
ain n is best taken as the largest integer less than 1/b; i.e., we choose the minimum number of 
al Servicing intervals. Note that in this case the expression for f(x) need not be known. 


Problem 3 
As an introduction to problems involving random deterioration functions, we will consi- 


no longet der certain preventive maintenance problems. In Problem 3 it is required to keep a certain 





rviced to System in operation for as long a time as possible without failure. Suppose the operation of the 
ned to be System depends upon whether a particular component is in good working order. Further, 
ce). Afte assume that in use this component deteriorates until it eventually fails. A preventive mainten- 
ance policy is under consideration which involves replacing the critical component every tg units 
of time until a failure occurs. 

Suppose that the cost of a system failure is given by (d-b7), where 7 is the time lapsed 


until the failure occurs, and the cost of replacing the component is denoted by "a." Let N 





ns''! 


a 


» ‘epresent the number of replacements before the system fails and L the length of time taken by 


os 
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a new component to become defective. (Assume that L has the distribution function Fr and 
density function f,)- 


The cost associated with a policy of replacing the component every t, units of time is 
C(t,) = a N(tp) + d-br(to), 
and the expected cost is 
(3) EC(t,) = a EN(ty) + d-bE7(to). 
The relationships between the distribution functions of 7, N and L, as well as their expecta- 


tions, were obtained by Weiss [ 25] . 


The density function of 7 is 


k 
f(t) = [1-Fy(to)] £1 (t-ktg), [ktp <t < (k+1)to], 
where k = 0,1,2,.... 
Hence, the expected time to failure is 
" (k+1)t, 
k 
Er(to) = >) \ t[1-F,(to)] dF, (t-kto), . ( to > 0.) 
k=o “kt, 


After manipulation, this reduces to 


’ 
[ 1-F, (t)] dt 
Ext.) 22 L 
F, (to) 





The probability function of N, obtained from its relationship to 7, is 
(k+1)to 
P }N=k{ -| f _(t)dt, k=0,1,2,... 
kt, 
Hence, the expected number of replacements is given by 
(k+ 1)t, 


EN(t,) = du f_(t)dt (tp > 0) 
~ 2 


kt, 


“dy k [ 1-F, (to)]* F (to) 


1-F} (to) 
Fy (to) 


The expected cost, from Eq. (3), is 


f° 
1-F (ty) ee i [1-F,(t)]at 


Fy (t,) Fr(t)o 





EC(t,) =a 


Unfortunately no simple characterization of Fy is apparent which will inform us when the 
value of t, which minimizes the expected cost is either finite or infinite. 
Suppose we now consider the same preventive maintenance problem, except that 


the replacement process is repeated indefinitely. We note in passing that an obvious 
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counterpart to this problem exists for items which deteriorate in storage. 
; Problem 4 

- To review: the replacement policy under consideration is as follows: 
If an item fails before t,, replace it. 

If the item does not fail by tp, replace it. 


We assume the cost structure 


ta- d'-bT; fH 0= Tj < to 


C( Tj» to) = 


to], 


where b> 0, d'=a+d, and7 j; is the length of time that the it? item is in operation. Again, if 


it is assumed that the random variable L,, denoting the time taken by a new item to become 
defective, has a known continuous distribution function Fi the expected values of 


0.) 7; and C; are, respectively, 


1\/ 


to 
ET=ET7, “| tf, (t)dt + to [1-Fy(t,)], wid... 
and 


EC 


to 
ECj -{ (a'-bt)f, (t)dt +at,[ 1-Fp(t,)], i=1,2,.. 
) 


Over an interval of length T, the total cost is 


N 
G(T, t,) “iC, 
1= 


and its expected value is , 
EG=E (2) 


i=1 


where N is a random variable representing the total number of items in use during the 
interval. 
It can be shown that the process under discussion is mathematically equivalent to a 


Sequential sampling process which terminates at the smallest value of N for which 


) 72T. This equivalence enables us to use a fundamental result of sequential analysis due 
i= 


to Wald [ 24] to compute the expected value of the total cost according to the formula. 
EG = EN-EC. 
Since the process is repeated indefinitely, we are interested in the average value of 
the expected cost per unit time, denoted by G, and given by 


EN:EC . 
T 


the 
G(t,) = lim 
T7370 
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A well-known result of renewal theory (see, e.g. | 23]) is that the expected number of 
























replacements for T very large is approximately 


Thus, we obtain Gt,) ” a? 
As in the previous example, we do not have simple criteria for determining when the cost mini- 
mizing value of ty is finite. 

Another, well-known, replacement problem was described by Campbell [ 4] . 

Problem 5 

A quantity (n) of new items are put in storage. Assume (i) that the items deteriorate 
while in storage and (ii) that the times taken for each item to become defective are independent, 
identically distributed, random variables with distribution function F;. As the items become 
defective they are replaced with new ones. After an interval (t,) all n items are replaced, 
regardless of their conditions. 


Suppose that the cost of replacing an item before tg is equal to b while the cost of 


replacing all at one time is equal to a. Let N(t,) denote the number of replacements (failures) ‘ 
up to time to. The problem is to find t,) so as to minimize the expected average cost per unit f 
of time. If C denotes this cost, it may be represented by the expression V 
I 

1 

C(to) = i [a + bN(to) | ls 1 


I 
Hence, the function of interest, its expection, is given by 
1 
EC(to) = t [a + bEN(to)] . I 
. t 
The computation of an optimal value of t, first requires solving for H(t) = (1/n) EN(t) ir 
in the renewal equation t 
fe) 
EH(to) = Fy (t,) +f. EH(to-t)dF (t). 
ss 
Asympotic (i.e., as to—) results of a fairly general nature are available, (See 
Smith | 23] .) However, for the problem in hand, we require knowledge of the transient 
behavior of EH(t,), which may be secured only by making specific assumptions on the form 
of Fj. For example, suppose F, has the form of a gamma distribution with parameter su 


a@=1and6 =1/u, i.-e., 


UX 4. (t>o) 


t 
F y(t) | u2xe™ 
Oo 
Bartlett | 1| has shown that 


1 a 
EH(t,) = 5 [uty - e *o sinh(ut,)] . 





lini- 


ate 
dent, 
me 
ed, 


ures) 


unit 


t>o) 
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Hence, -ut 


_a , bn . = 
EC(t,) = to +> E - S" Samat) 
The unique value of tp which minimizes EC(t,) satisfies the equation 


(1 + 2ut,) -(1- 2a )_2uto 


provided that a<bn/2; otherwise, the expected cost is minimized by the policy of replacing 


defectives as they occur, i. e., tp =~. 


INVENTORY DETERIORATION 

In this section we shall describe twogroups of problems which occur when items 
deteriorate in storage. As contrasted with the usual inventory problems involving the deter- 
mination of optimal stock ordering policies, the first group is directed towards finding optimal 
stock issuing policies. The second problem area deals with the development of inspection 
procedures to maintain a prescribed level of quality in storage for these deteriorating goods. 
Stock Issuing Problems 

Problem “ 

At the start of a certain period there are N items of the same kind in storage; their 
ages are denoted by S;,...,Sn- The deterioration of each item is represented by a utility 


function Uj(t), i=1,...,N. The demand for the items is reflected in a schedule (a,, aioveeaall 
th 


) 
which exhausts the stock, where aj denotes the time from the present at which the i nity 
must be issued. It is of interest to determine the issuing policy, i. e., the order in which the 
items are to be issued, which maximizes the total expected utility from the inventory while 
meeting the specified schedule of demands. 

Case I.° Suppose that the analytical form of U;(t) is known. It is easy to show that this 
issuing problem can be formulated as a special kind of linear programming problem known as 
the "assignment problem." Let ujj denote the utility resulting from the issue of an item of 


initial age S; to meet the demand at time aj, i.e., 
Ujj = Uj(S; + aj) G;. j= 2,0... 


Then the problem requires the maximization of the linear form 


N N 
> > UjjXij; 
j=1 i=1 
Subject to the restrictions 
xij = 0, 1, hj 3% dys 


N 
= * z, Xj = d. 


—_ i=1 j=l 


Me 


, Greenwood [12 ] should be credited for calling attention to this problem area. 
Discussed in [8] by Derman and Klein. 
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Computational procedures suitable for the solution of the assignment problem are given by 


Kuhn [ 14] and Munkres [18]. It may be noted that the same computational problem exists 






either when the items have different ages and the same deterioration function or when they 
deteriorate at different rates and the ages are the same. The same formulation also holds 
if the stock ordering schedule is known for a number of additional items and the depletion 
schedule is extended to include this additional quantity. It is interesting that in some situa- 
tions (described below) it is not necessary to know the analytical form of the deterioration 
function. 

Case II. Suppose it is known (i) that U,(t) = U(t), (ii) that U(t) is a convex function and 
(iii) that S} << Sp< ...<S,. Then it has been demonstrated [ 7] that an optimal stock issue 
policy is one known as LIFO (last in, first out), i.e., the youngest item is issued first. Note 
that only knowledge concerning the relative ages of the items is required for this case. 
Further, the withdrawal schedule need not be known either. If U(t) is a concave function the 
opposite issuing policy — FIFO (first in, first out), i.e., the oldest item is issued first — is 
optimal. 


Case II. Suppose that Uj(t) = U(t), S1<S2< ..<Sp and that items are to be continually 
added to stock and continually issued. Then it has also been shown | 7] that 


(a) If U(t) is convex, the best issuing policy is LIFO, and 
(b) If U(t) is concave, the best issuing policy is FIFO. 
Again we note that the withdrawal schedule need not be known for these results to apply. 
Problem 2 
There are circumstances in which the demand schedule will not be independent of the 
deterioration function as it was in the above. Suppose that the deterioration of an item in 
storage is reflected in the length of its useful life when it is issued. When the item’s field 
life is over, a replacement is issued from inventory. Assume that we are given the relative 
ages of the items at the start of this process and that all of the items deteriorate according 
to the same law. Since LIFO and FIFO policies are easily understood and implemented, it is 
of interest to know the conditions under which they are optimal. 
Let L(S) denote the field life of an item which is of age S when it is issued. Then, the 
following results are known (See [7] and [15]): 


THEOREM Ia: LIFO is an optimal issue policy for a stockpile of n>2 items if either 
of the following conditions hold: 


(a) L(S) is a convex function and LIFO is an optimal policy for n = 2. 
(b) L'(S) >- 1 and LIFO is optimal for n = 2. 


Examples of deterioration functions for which LIFO is an optimal issue policy are: 





(a) L(s)= >, (a>0, b>0) 





(b) L(S) = ce-kS (c >0, k >0) 
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THEOREM Ib: FIFO is an optimal issue policy for a stockpile of n> 2 items if 
either of the following conditions hold: 
(a) L(S) is convex and FIFO is the best policy for n = 2. 
(b) L'(S) >-1 and FIFO is optimal for n= 2. 
THEOREM 2: _ FIFO is an optimal policy if 
(a) LS) >- 1, 


(b) L(S) is monotone 9 and concave. 


The characterization of deterioration functions for which LIFO and FIFO issuing 
policies are optimal, as given in Theorems 1 and 2, is essentially unsatisfactory. 
Although Theorem I implies that the choice between LIFO and FIFO may be made largely on 
the basis of a solution to the two-item problem, the subclass of deterioration functions for 
which either of the policies is optimal for this case is unknown. Similarly, despite the fact 
that Theorem 2 gives conditions which do not require special consideration for n = 2, 
condition (a) probably can be weakened since if L(S) is linear (with any slope) FIFO is still an 
optimal policy. Further, necessary conditions for optimal policies to be either LIFO or FIFO 


are also unknown. 


SURVEILLANCE SAMPLING: 

Derman and Solomon [9] have considered the problem of maintaining a preassigned 
quantity and quality level in stored items which deteriorate. To accomplish this end, they 
suggested the use of rules involving periodic sampling of the inventory. 

Consider one lot of items. Suppose this lot enters storage at time t=0 and itis of 
quality p(0), i.e., the initial proportion defective is p(0). Further, assume that during its 
stay in inventory it deteriorates according to some known law. Then, by the use of a decision 
rule which at any time either permits a lot to remain in storage or, by rejecting a lot, 
simultaneously calls for the entrance of a new lot of quality p(0), the quality of the inventory 
can be kept (in probabilistic terms) above a preassigned level. In general, the details of 
such a rule (e.g., the severity and the frequency of inspection) may be made a function of the 
age of a lot to obtain increased protection against undesirable eventualities. 

Using this approach, the age of the lot in storage at any time t can be viewed as a 
Markov process with a countable number of states and with stationary transition probabilities: 

Pjyj+1 = probability of a lot going from age j toage j + 1; 


= probability of accepting lot of age j + 1 onthe j + ith 
inspection; 


probability of discarding alot of age j for anew one on the j + yth inspection; 


Pj, 0 
1 Pj, jth 
If now, instead of one lot, we assume that there are a large number of lots being 


handled simultaneously and further that the inspection-replacement process has been operative 








Prof. G. Lieberman has informed us that P. Zenna, Stanford University, has shown that 
the monotonicity requirement may be dropped. 
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for a long time, then the Markovian steady-state conditions together with the law of large 
numbers may be assumed to be applicable. Thus, the expected proportion defectives 7(t) in 


the system at any time t between inspection intervals is given by the expression. 





(t)= > pitty, (o< t <1) 
j=0 


where p(j + t) = proportion defective in a lot at time t after the jth inspection (obtained 


from knowledge of the deterioration function) 


and Vj steady-state probability of finding a lot of age j in the inventory system 


The values of vj are obtained by solving the system of equations 


Vj+1= VjPj> f+ b (j=0, 1, ----) 
= vj = a 
j=o 
subject to the constraints 
Vi > 0. (j=0, 1, ---) 
a3 
They are J 
vj = Yo bid Lr + i[p(r + 1)| 
r=0 
d 
an - -— ¥ 
Vo=J1+ DU Lr +ilpr+1)} , 
j=1 r=o 
where Lj | p(i)] = probability of accepting a lot of age i if the proportion defective in the 


lot is p(i). 
This latter quantity is the operating characteristic of the surveillance sampling plan. 

By virtue of the above expressions, it is possible to evaluate the performance of any 
surveillance sampling plan in terms of 7(t); hence, plans can be chosen in ways analogous 
to those used in acceptance sampling. Detailed work is required in this area. 

More recently, Frank [11] considered a continuous analog of the same situation. 

He allowed the interval between lot inspections to be continuous random variable with a 
known distribution function instead of a constant. The main effect of this change from the 
discrete model just discussed is that the expected proportion defective is a continuous 


function of time and in the steady-state it becomes a constant. 
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The purpose of this paper is to illustrate the applicability of two of the 
fundamental tools of classical analysis, functional equations and succes- 
sive approximations, to various classes of multidimensional maximiza- 








tion problems. 


The advent of the modern digital computer makes it particularly appro- 
priate that this synthesis of new and old be made. Today we can envis- 
age as single steps in an iterative solution the resolution of particular 
problems which even tenyears ago would have beenconsidered as beyond 
the scope of mathematical power. Since the rules of the game have 
changed so drastically insofar as what we accept as a feasible solution, 
it is essential that we revise many of our techniques with this well in 
mind. 


Suitably combining classical and modern devices, we can hope to make 
a start on the solution of the many complex and significant problems that 
challenge us. 











INTRODUCTION : 
The purpose of this paper is to illustrate the applicability to two of the fundamental 


tools of classical analysis, functional equations and successive approximations, to various 





classes of multidimensional maximization problems. 

The advent of the modern digital computer makes it particularly appropriate that this 
synthesis of new and old be made. Today, we can envisage as single steps in an iterative 
solution the resolution of particular problems which even ten years ago would have been con- 
sidred as beyond the scope of mathematical power. Since the rules of the game have changed 
so drastically insofar as what we accept as a feasible solution, it is essential that we revise 


many of our techniques with this well in mind. 
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Suitably combining classical and modern devices, we can hope to make a start on the 
solution of the many complex and significant problems that challenge us. In this paper, we 
outline the mathematical ideas. Computational results may be found in a number of the refer- 
ences and in a forthcoming book by S. Dreyfus and the author. See, in particular, [3], [4], 
[8], [9], [12], [19], [20]. 


A FUNDAMENTAL PROBLEM 
The various problems upon which we shall focus our attention can all be regarded as 


special cases of the following question: 'Determine the maximum of the separable function 
(1) F(x 1,Xg, ---»Xy) = (X41) + Golxg) + --- + EnlXy), 
subject to the constraints 


N 
(2) 2 aii) <cj, i=1,2,...,M, 


and, usually, the nonnegativity condition 
(3) oO, 5089. ...8* 


Naturally, we don't expect a uniform solution to a problem of this generality. 


DISCUSSION 

The problem is usually trivial from the standpoint of the existence of a solution. In 
most cases, the x; are subject to constraints which allow only a finite set of choices. Con- 
sequently, the problem that we face is that of using mathematical ingenuity to determine the 
maximum value by means of something better than an enumeration of possibilities. 

It is, first of all, desirable to obtain as many general methods as possible which will 
yield a solution, by means of hand computation or machine computation, regardless of the 
form of the functions, gj(x), or the constraint functions, ajj(x). It is the great merit of the 
simplex method, and recent improvements, that in the special case where the functions 
g,(x) and aj j(x) are linear, where the number of variables and constraints are not too large 
and where continuous variation is permitted, it serves this useful purpose. 

Similarly, if there is a "small" number of constraints, dynamic programming furnishes 
a solution, regardless of the form of the functions and the type of variation that is permitted. 
These, however, cover only a small part of the domain of problems. It must be emphasized, 
then, that at the present time "universal solvents" do not exist, and it is highly unlikely that 
they ever will. Consequently, various special devices must be employed in which both the 
analytic structure of the equations and the intrinsic structure of the physical process giving 
rise to these equations are utilized. 

Amusingly enough, there are some who feel that this is a form of "cheating."" There 
is an interesting discussion of this point in N. Wiener's autobiography, in which he comments 


on the difference between Paley and himself in their attack of a mathematical problem. 
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This idea of divorcing mathematical research from the physical origin of the problem 
runs counter to the whole stream of scientific development. It not only leads to decadent 
mathematics, it leads to lame and limping analysis. 

It should be realized that the mathematical equations that one obtains in discussing a 
physical process are not equivalent to the physical process but only one projection of the pro- 
cess, in a special coordinate system, and from a special point of view. To disregard the orig- 
inal process once the equations have been derived is to throw away valuable information. 

It is more than likely that when mathematics has reached a sufficiently high level of 
sophistication the solution to many of the basic problems which today are treated in complex 
analytic terms will be given in simple verbal terms. An approach such as that contained in 
Boldyreff's "flooding technique" is a presage of developments of this nature. 

All this is not to say that analysis of the type that follows is superfluous. Only by 
attacking and solving a variety of individual problems can we hope to discern the general pat- 


terns to the solutions that we hope exist. 


DYNAMIC PROGRAMMING APPROACH - LINEAR FUNCTIONS 


We shall begin by considering the problem of maximiziang the linear function 
(4) L{x},Xo; ints » Xp) = byx, + boxg + .-- + Dyxy, 


subject to the linear constraints 


(5) Di S¢ coe | 


and the nonnegativity requirement x; > 0. The coefficients ajj are taken to be nonnegative, 
with a sufficient number positive so that the x; are uniformly bounded. 


The maximum value is a function of the Cj and N. Hence, we write 


(6) fyylC Co, - +++ Cyg) = Max L(xy, XQ,---Xy)- 
xX; 
i 


For the case N = 1, we have the result 
(7) f4(Cy,Cg---»Cyq) = Max byxy, 


where the maximum is taken over the region 





(8) 24%, < cj, i=1,2,...,M. 
Hence, 
—(Pici Pace bicyy 
®) £4(Cy,CQ,---,Cy)= Min| -—— , 3 ge ss 
411 421 ami 


Using the principle of optimality [1] for N > 2, we derive the recurrence relation 


(10) fle 


12Cg>+-+9Cyq) = Max [Dy Xp +fy-y 1 - 8yyXyoee uM -4@qn*n)]° 
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where xy is subject to the conditions 


(11) O <x, < Min 


This procedure furnishes an iterative algorithm for determining the solution of the 
original maximization problem. Starting withthe known function fj (Cy»Cos Heng Cry)» 


we determine f then fg, and so on. In so doing, it furnishes the solution to a class 


, 
of problems, ‘omiietins to different values of the ci and N, Thus, automatically, 
this method yields a stability or sensitivity analysis. 

Computationally, we face grave obstacles when M is large. For M = 1, 2, or 
even 3, or in cases where the c, can assume only a small set of values, we have a 
feasible algorithm performed on high-speed electronic computers. For large values of M, 
we face the impossible task of tabulating functions of M variables. For these values 
of M, we must admit that this method in its lowest-level most routine form is not workable. 
In what follows, we shall discuss various techniques that can be used to remedy some of the 


difficulties. Fortunately, the most interesting problems lie ahead. 


CONSTRAINTS 

When applicable, the algorithms of linear programming, such as the simplex method, 
are very powerful. Problems of large dimension can be rapidly solved. In general, however, 
the use of these techniques is predicated upon the assumption of continuous variation of the 
variables Xs in their regions. 


Integral constraints of the form 


(12) x, = 0, 1, 2, oss, 


of which the extreme form is x, = 0 or 1, _ introduce grave complications. Furthermore, 


the presence of additional constraints, such as 


(13) <a < & 
Ix i — 1 
or 
(14) aS Xiy1 -%1 LD, 


can seriously impede the rapid determination of the solution. 


"Complementary" constraints, such as 


i = 1,2,---, N—1, 
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again escape the routine application of the standard algorithms of linear programming; see 2]. 

Restrictions of this type, however, greatly increase the efficiency of the dynamic pro- 
gramming technique - when it is applicable at all. Generally, the more restrictions on the 
freedom of variation of the xj, the more rapid the dynamic programming solution (See the 
exercises in [1] for many examples). 

Finally, one knows, in advance, that computing time required to yield the solution to 
the original problem, plus the problem for a range of C values and values of n in the range 
1, 2, ---, N, is directly proportional to N. In the use of linear programming the expected 


time for obtaining a solution is of the order of magnitude of Nn? or n°, 


REDUCTION IN DIMENSIONALITY 


Having pointed out how drastically the dynamic programming is curtailed by dimension- 
ality, it is important that we point out how the linearity (really the homogeneity) of the prob- 
lem can be exploited to reduce the dimensionality of the sequence of functions by one. This is 
quite significant, since it means that processes involving two constraints can be treated by 
means of sequences of functions of one variable. Combined with the Lagrange multiplier for- 


malism discussed subsequently, it results in a great increase in the scope of our methods. 


Consider, to illustrate the idea, the problem of maximizing the linear function 


(16) U(X), Xo ++*> Xn) = b)x, + boxy + «++ + by ¥n > 
subject to the constraints 

N 
(17) be j * ", < Cy: 

N 

1 a= °° 


and x > 0. 


As we know, this leads to the functional equation 


18 = ~ np 

(18) fryley Cp) — [Pn * Ta « 4in*y? 2 n*n)) , 
N 

where xy is subject to the restrictions 


(19 Ci; C2 
< xy < Min|s—, aco 
) 0 < xy < Min 21N’ 20N 
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we 
It is clear from the original variational problem that each of the functions fy(cj, cg) must be 
a homogeneous function of c, and C9. Hence, for C1, Co 0, (o 
“9 ‘y 
(20) fy (Cy >a) = Cif i? c ,= Cofy oo ; he 

pu 

Consequently, setting 
Co NC 
(21) Xy = 2Cy, 6 = ey r 

val 
we see that (18) becomes 

(2¢ 

Cc. = aor 
e 2 2N 
(22) fy (7) = Max byZ + (1 — a,y2) fy-1 (, =~ ) 

Zz 1N wh 
where z_ is subject to the restriction (2¢ 
(23) 0< 2 < Min(2— a). 

1N 2N 

Pr 

Setting - 

tio 
(24) u(r) = f,(1,r), O< r < @ fu 

N N > , a ps ’ 

mi 
we see that (22) yields a recurrence relation for the function Uy (r), a function of one ca 
variable. 

It does so, however, at a certain expense. In the original two-dimensional formulation, we 
if O< cy < dy » OK Co < dy initially, then at no time do we have a point outside this 
rectangle. As a matter of fact, the grid is actually decreasing insize as N decreases. 7 

The one-dimensional formulation requires the tabulation of the functions Uy) over | 
some subset of the infinite interval 0< r < ~. This is a quite unpleasant feature, with fel 

aaa e 
all kinds of subtle opportunitites for error. We can get around all of this by computing the two | 
no 
4 \ ) elt eee 
sequences, } fy Gr); >) fy tr, 1); » both in the fixed range 0< r << 1. the 

Thus, we use (22) when z is such that 

pr 

a x, _ 4 
(25) Co Bont <1- Zain 

(3( 
When z is outside this interval, i.e., 

su 


(26) Cc. — &.F > = Zay ny? 
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this 
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1 — Zain 
= oe —______. ] 
(27) fy r) — byz + (Co aont)fy_y ip tae 

For an application of this idea to the study of "bottleneck processes" insofar as com- 


putational aspects are concerned, see [3]; for the analytic aspects, see [1]. 


NONLINEAR FUNCTIONS 
The same reasoning used in the discussion of the linear case shows that the general 


variational problem, posed at the beginning of this paper, leads to the functional equation 


(28) fylCy,Cgs--+5¢yg) = Max fbyxy + fy [ey — ayy ly) +--+ Cy — avnley)] } > 
*y 
where xy is now subject to the constraints 


(29) (a) Any) < Cy» gt ee eS 


(b) xy > 0. 


Provided that M = 1 or 2, or that the set of allowable c values for large M is not too great, we 
possess a Straightforward algorithm for obtaining the computational solution of the maximiza- 
tion problem. This method requires no information concerning the structural properties of the 
functions that occur, it allows variation over weird sets, and yet it yields the absolute maxi- 
mum. Consequently, it furnishes a basis for the computational solution of problems in the 
calculus of variations; cf. [4, 5, 6]. 

Withal, the question remains: ''What can be done when M is large?" To answer this, 


we shall dip into the storehouse of classical trickery. 


LAGRANGE MULTIPLIERS 


The first device we wish to discuss is the Lagrange multiplier. Although it is commonly 





felt that this technique is intimately associated with continuous variational problems, this is 
not the case. Actually, the many applications of this device are all reflections of moment 
theory. Since any amplification of this remark would take us too far afield, we shall merely 
present the formal treatment. 

Consider the problem of maximizing the function 
(30) F(x,» - 


g) 1219 %y) = 8y(X) + BolXp) + «+ + Bryn), 


Subject to the three conditions 
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N 
< 
2a, {(*,) < Cy 


N 
Pa 
2ag(x;) ae 


and the condition x; > 0. 
As pointed out above, this problem can be treated by means of functions of three vari- 
ables. Since this is one variable too many for current computers, we proceed as follows. In 


place of the foregoing problem, consider the problem of maximizing the new function 
N 
(32) F (x1) Xp) +++ Xpy) = 8y(Xy) + SolXq) + --- + SpylXy) — ‘2 9(xj); 


subject to the constraints 


Mz 


(33) Ail) < ¢, 
N 
< 
and x, > @. 


The quantity 4 , assumed fixed for the moment, is the Lagrange multiplier. The new 
problem leads to the sequence of recurrence relations 
\ 
(34) Oyy(Cy» CQ) = Max ) Sy(Xy) — Aayn(Xy) + Oy] [cy — ayn(Xy), Cg — agn(xy) | , 
x 
N 


where Xn is subject to the constraints 


(35) 21 (Sy) < Cy» Aon (Xn) < C9; xy 2 0. 


It is clear that A occupies the role of a price, an exchange parameter. Consequently, 


we vary A from 0 to ~until the third constraint is met. 


DISCUSSION 

There are a number of interesting aspects to this procedure. In the first place, ob- 
serve that we have not stated explicitly that continuous variation is required, and it need not be. 
As stated above, we are in reality exploiting a simple aspect of moment theory, which in the 
continuous case reflects itself in a duality between point and line coordinates; see [7] and [1] 
for further discussion and applications. 

Second, although the rigorous details connected with this method can become thorny, it 


is quite easy to prove the following result: 
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Let x0), X_(A); breton xy be a solution obtained from (34) and compute the quantity 














(36) N 
DX a,.(x,(a)) = ¢,(a). 
—-« ob & 3 
i=1 
Then [x,0), X_(A); ee » Xa (A) ] is a set of values yielding the maximum of the function in (30), 


subject to the constraints 


(37) 


Mz 


a, (%) Se, 


vari- is] 
3. In 
N 
& Ag,(x;) < Co, 
N 
2 ag (xi) < Ca(A)- 

It follows that as we vary A between 0 and ~, we obtain a set of Cy values which we hope 
cover the range of interest to us. It is clear that as A increases, the value of Ca(A) decreases 
monotonically. Consequently, one can employ fairly efficient search techniques to ascertain 
the appropriate range of values for A. Some applications may be found in [8], [9]; the 
original discussion is contained in [10]. 

The merit of the method lies in the fact that it enables us to trade time for memory 

e new capacity. If we have quite limited capacity, we can introduce two Lagrange multipliers, and 
consider the problem of maximizing the function 
. N N N 
(38) Fo(XysXpo +++ Xy) = 2s g(x) —A De agi(x;) — AD agi(x,), 
i=1 i=t . or 
subject to the constraint 
N 
(39) 2» a, ,(%) < Cy x, >0. 
iently, 
Furthermore, given no memory capacity, we can attack the problem of maximizing the 
function 
N N N N 
sie (40) F(x) X,, wins »Xy) = ~ g,(x,) —>Ay 2 a, (x) — A» Pm Ao,(x;) — rz 2 ag;(x,). 
1 not be. The point is that by introducing parameters from the dual space, the space of the Ai 
in the and working partly in one space and partly in the other, we have a means of decomposing an 
nd [1] M-dimensional problem, one in which M constraints occur, into a sequence of (M — k)- 


dimensional problems. The arithmetic of computers is such that this is quite effective. 
orny, it 
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SUCCESSIVE APPROXIMATIONS 

The method we have just discussed may be considered to be an application of the 
technique of successive approximations, since we guess a set of Ay values, solve the resulting 
maximization problem, test to see if a solution has been obtained, and continue in this fashion, 

In what follows, we wish to present some direct applications of this basic technique of 
analysis. Before so doing, let us make some preliminary remarks. There are really two 
problems to be considered in dealing with physical processes which lead to mathematical 
problems of the type we have been considering. The first is that of determining absolute 
maxima or minima, the second is that of obtaining a feasible policy which is significantly 
better than current policy. 

It may frequently happen that if we find a feasible policy of this type quickly, we may do 
more good than with an optimal policy found too late. 

In the musical comedy, "Annie, Get Your Gun," there is a charming duet between Annie 
Oakley and Frank Butler to the theme of "I can do anything better than you can."' Our aim is to 
present various methods which enable us to start with any given policy and improve it succes- 
sively. These steps may or may not lead to the actual optimum, but they do yield monotone 
improvement. We shall not discuss the difficult convergence questions here. It is to be ex- 
pected, as in the application of the usual gradient techniques, that local maxima and saddle- 


points of more complex type will cause a certain amount of grief. 
SUCCESSIVE APPROXIMATIONS AND LINEAR FUNCTIONS 
Let us return to the problem of maximizing the linear function 


(41) L(x) = L(x), Xp, gis » Xp) = bX + DoXo i oes Dyn? 


subject to the linear constraints 


N 
(42) 2 ajjXj S Cp» = £26055 
j=1 
and Xj Pa 0. 
0 0.0 0 ie aa , 

Let x = X9Xgo++ +> Xx be a point inside the polyhedral region. To start our 
approximation, let grees xy be fixed at the value i worry xy and consider the problem of 
maximizing 
(43) LX)» Xp) = byX; + boXp, 
subject to the constraints 

™ 0 
(44) AjyX] + AjgXq < Cj — = i= 1,2,...,M, 


and x; > 0. 
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This problem can be tackled in several ways. First of all, we can employ a straight- 
forward search technique, fixing a value of Xy> maximizing over the allowable set of Xo values, 
(which is trivial), choosing another vaiue, and so on. Second, we can employ the dual formula- 


tion, "Minimize 


(45) M 
ie CiYi> 
i=1 


over all y; satisfying the constraints M 


M 
(46) > ai1yi > 1 
i=1 
M 
i=l 
Third, one can use the simplex method itself. 
Let 
(47) xi = x! xd x0 0 
1’ 2? 3? eeey *y 


be a set of X; values obtained in this way. It is clear that 


0 1 
(48) L(x’) < L(x ). 
Now fix the variables x,,x X,, at x! x? ° and consider the problem of 
} 17% 4707+ xy o arr xy 
maximizing boXo + byX3 over the allowable set of (Xo; X3) values. Call the resulting set of 


values x . Since x, =X,, X, = X, is an allowable choice of x, and x,, it is clear that 
1 9 2 2 3 3 2 3 
L(x )< L(x). 
Continuing in this way, we obtain a monotone increasing sequence of values of L(x). 
Provided that we started sufficiently close to a set yielding the maximum, it is clear that we 
must converge to it. One would suspect that the absence of relative maxima in this case will 


yield convergence in general. 


STOCHASTIC CYCLING 


As outlined inthe previous section, we have envisaged a steady progression, first 


x,» Xo), then (Xo, Xa), then (xg; X4), and so on. In general, to avoid individual idiosyncrasies 


of problems, it may be better to proceed in a random fashion. Let the digital computer 
generate a random integer between 1 and N at each stage, say i, and use the pair (x; xy 1) 
It might also be worth occasionally, in a stochastic sense, varying over a triple of 


values (xj>X, | 1 *i+ 9) or over a double random set (x,, x;)- A great deal of mathematical 


experimentation of this type remains to be done. 
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SUCCESSIVE APPROXIMATIONS AND NONLINEAR FUNCTIONS 


It is clear that precisely the same ideas can be used to treat the general nonlinear 
maximization problem given at the beginning of this paper. Since, in general, there will be no 


dual problem, a straightforward search technique will probably be best. 


HITCHCOCK-KOOPMANS TRANSPORTATION PROBLEM 

Quite often, in place of applying successive approximations in this uniform way, inde- 
pendent of the particular problem under discussion, it is far better to use the structure of the 
underlying process as a guide. 

In the Hitchcock-Koopmans transportation process with nonlinear cost-functions, we 


face the problem of minimizing a function of the type 


(49) F(x; ;) Py di (3 5)» 


subject to the constraints 


(50) Bay Sy coe ee 


i vj» j= 1,2,...,N, 


and x..>0, where 2s; = =r 
| eae i j 
Here Xij represents the quantity sent from the "source" i to the "sink" j, and di, (x, ) 


% 


represents the cost of transporting this quantity. This function, if taken to be linear, leads to 
a problem for which many elegant solutions exist. If, however, di, (2) is a nonlinear function of 
z, many of these techniques fail. 

To treat the problem by means of successive approximations, we fix the allocations from 
M — 2 of the sources and minimize over the allocations from the remaining sources. It turns 
out that by virtue of the condition Ys. = Zr. this leads to a problem involving sequences of 
functions of one variable, rather than > 5 a of functions of two variables. The details may 
be found in| 11]. 

Computational experimentation with this method has led to satisfactory results [11]. 
What we wish to stress is that the physical process underlying the mathematical problem leads 
us to a very natural type of successive approximations, a type of approximation which the 


mathematical problem by itself might not suggest. 


BROWN-VON NEUMANN FICTITIOUS PLAY 

Another very powerful technique hinges upon the idea of regarding the maximization 
problem as the limiting steady-state version of a multistage maximization process. This idea 
was applied to the determination of the value of games by Brown and von Neumann, with a 


rigorous proof first supplied by J. Robinson [12]. The method has been brought to a much 
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This concept is also the essence of the gradient method; see the article by Rosenbloom [13] . 


QUADRATIC PROGRAMMING 
As an application of successive approximations in an entirely different way, consider 


the problem of maximizing the quadratic form 


N 
(51) (Bx) = a bir 


subject to the set of constraints 
N 

(52) > aijX < Cy 1=1,2,-.-,M, 
j=1 


and x; > O; £= B523«+05- 
Let us discuss the case where B is positive definite, which means that the maximum 
will occur on the boundary of the region described by the foregoing constraints. 


Let us start with observation that 
(53) (x.Bx)M|y+x—y, Bly+x-—y] 
= (y, By) + 2[(x — y), By] + [x—y, B(x —y)], 


which means that 
(54) (x, Bx) > (y, By) + 2(x — y, By), 
for all y, with equality only if x= y. Hence, 


(55) (x, Bx) = Max [(y, By) + 2(x — y, By)]. 
y 


Consider, in place of the original nonlinear problem, the problem of maximizing the 


linear function 


(56) Pa tte ~ 


over all x, subject to the constraint of (52), where x? is an initial approximation, satisfying the 
constraints of (52). 


Let x! bea maximizing set, xl=(x1, x}, eee Xf) . Then, from the maximizing property, 


(57) (2° px) + a(x! - i, Bx”) > (.°, px) + (x9 - x?, Bx”) 2 (x°, px?) ; 
Furthermore, from (54), 
(58) (1 ; Bx?) > (0 ‘ px’) + a(x! “- ; Bx°) ; 


Combining these last two results, we have the inequality 


(1 < px!) = (0 ’ px?) ° 
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Repeating this procedure, with x! taking the place of x’, we obtain a set x? yielding the in- origin 





equality of hig! 
multis 
2 2 1 1 
(60) (x", Bx”) > (x’, Bx) > (x?, Bx), wes 
these | 
and soon. We therefore have the desired monotonicity. The question of convergence is a 
more difficult matter which we shall not discuss here. lat 
relatic 
CONVEX FUNCTIONS IN GENERAL (66) 
The foregoing technique can be applied to a number of other problems, since the basic 
: yhere 
property we are exploiting is convexity. If g(z) is a convex function of z for all z, we can 7 
write (67) 
(61) g(z) = _ | gly) + (z - y)g'(y)]- 
Proceeding formally, we see that the function g(y) + (z — y)g'(y) has a stationary point at the —* 
value of y for which the derivative vanishes, ; functic 
(62) (z — y)g"(y) = 0. (68) 
If g"(y) > 0, the only solution is z= y. The result can, however, be established under milder 
assumptions. provid 
Consequently, if the functions gj (z) are all convex, we have a simple way of obtaining a discus 
sequence of monotone approximations to the solution of the problem of maximizing by one 
(63) F(Ky,X9,-++)Xyy) = 84 (Ky) + Soo) + --- + Ey Ky), — 
subject to the constraints mation 
N 
(64) 2 ae S Cj. i=1,2,...,M, ANAL’ 
and x; 2 0. (69) 
Similarly, if the function of N variables, F(x;,X9, ice Xn), can be written in the from 
N oF n= 0,1 
(65) F(x,, Xor+++ 9X) = Max [Py ¥py-+2y) # 2X _ y;) 4 ’ 

y i=1 (70) 
the foregoing techniques can be employed. Functions of this type are called Schur-convex Then 
(cf. Ostrowski |14]) and play a role in the theory of inequalities where such representations 
are extremely useful. Cf. Beckenbach-Bellman [7]. (71) 

whence 
MULTISTAGE PROCESSES 
(72) 





So far we have been using the functional equation technique of dynamic programming, 





n 
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originally devised to treat multistage processes, as a device for converting static processes 


of high dimension into dynamic processes of low dimension. In many physical situations, 





multistage processes arise naturally. This is particularly the case in the study of production 
processes, Of which "bottleneck processes" are an interesting type. A detailed discussion of 
these processes may be found in | 1]. 

The analytic problem reduces to the following. "Given a vector-matrix recurrence 


relation of the form 
(66) Xn+1= AnX¥n + Yn» X0=¢; 
where the y, are constrained by conditions such as 
(67) (a) Biyy < CoXp 
(b) yy, 2 9, 


how do we choose the y,, and perhpas the A,, B,, and C,, so as to maximize the terminal 


function 
(68) J(y) = (xy, a)?" 


This problem may be treated analytically (see [1], [15], and, computationally, [3]) 
provided that the dimension of the vector c is not too large. As in the other linear processes 
discussed in the beginning of the paper, the dimension of the process can always be reduced 
by one by using the homogeneity of the equations. 

We wish to use functional equations in a different way here, a way which will permit us 
to exploit the structure of the process and which will enable us to apply successive approxi- 


mations to more complex processes. 


ANALYTIC PRELIMINARIES 


Let us begin by studying the linear inhomogeneous difference equation 


(69) Xn+1 = AXy + Yn» XQ C; 
n=0,1,.... To obtain a representation for X, in terms of the y,, write 
(70) =A"z 


Then, substituting in (69), we have 


‘ Aa (n+ 1) 


(71 = . 
) 2n+1 2 Yn? 


whence 


n—1 


(72) zzc+ DA “ty, n>. 


n 


k=0 








































78 RICHARD BELLMAN 


This yields the desired relation 


n =) n—k—1 
(73) n odes Bee Vic? n?>l. 
k=0 
Observe that n — k — 1 - 0 fork =0,1,...,n — 1, so that no inverse of A need be taken — 


an important point. 
It follows that 


(74) 


N-1 
Xn = b+ 2B 


where b and B, are independent of the y, and fixed as long as N is fixed. 


A MAXIMIZATION PROBLEM 


Consider now the problem of determining the y,, subject at the moment only to con- 


straints of the form 


—1 
(75) (a) pan <k,, 
k= 


(b) ay S Vix S By 
(where yj, is the jth component of y;,) so as to maximize a prescribed function 


(76) g(Xjn> XOn> eoey XkN) 


of the components of xy. 

In mathematical economics, we encounter maximization problems in connection with 
multistage production processes and minimization problems in connection with scheduling and 
smoothing processes. In engineering control processes, we are usually trying to minimize 
the deviation of the system from some desired state. 


Introducing a Lagrange multiplier, we pose the problem of maximizing the function 
(77) 8(X4n> Xone +> Xin) -A ny 


We now wish to take advantage of the fact that in many situations of interest the number 
of components appearing in (76) is quite small. In any case, as we shall see, this type of pro- 
cess can serve, via successive approximations, as a stepping stone for the treatment of the 
more complex processes. 

Using the representation for xy in terms of the y, given in (74), we obtain the problem 


of maximizing 


N-1 N—1 = 
(78) g ( + Dobie --bp + Dy breve) — a hy) 
k=0 k=0 k=0 
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over the allowable vectors y,. We now regard this maximum value as a function of the 
scalar parameters bj, bg,..., by and N. 
Let us then introduce the new sequence of functions 


N-1 
(79) fyy(by, bg, ---, by) = Max a, -ee,Dp t+ ---) mA suo 
y k=0 


For N= 1, we have 


(80) {,(b,,bp, : -+yb,) = Max [a + bi 9Yo v+eyb, + by n—1Y9) _ sno) . 
yo 
while for N > 1, we have 


(81) fy(by,b9, eee , by) = Max - Ah(yg) + fy—1(b; + bi 0Yo> vee » D, + oor 
Yo 
Provided that r does not exceed 2, we have a simple, straightforward, computational 
algorithm. If r > 3, we must either restrict the range, use successive approximations as 


discussed below, or employ a still different method to be discussed below. 


APPLICATIONS 
This method has been applied to the study of the caterer problem | 16] (where explicit 
analytic solutions can be obtained in some cases) and can also be used to study inventory 


processes with time lags [17] where stochastic effects enter. 


QUADRATIC CRITERIA 

If the function h(z) is quadratic in the components of z, and, similarly, if 
AXi nN» Xn Biot Xn) is a quadratic function of the x;y, then the problem is one that can be 
attacked in a straightforward way by using the classical methods of calculus — provided that 
we impose no further constraints on the Yj- 

The temptation to do this, without further thought, must, nonetheless, be resisted! 

The variational equations will be linear, it is true. However, if the dimension of each vector 
is M, an N-stage process will give rise to a linear system of dimension MN. For M = 10, 

N= 10, this is the quantity 100, formidable but not overpowering; for M = 10, N = 100, this is 
overwhelming. 

By returning to (81), it is easy to show that fy(b, ...,b,) is a quadratic form in the 
bj, under the foregoing assumptions. Hence, (81) can be used to obtain recurrence relations 
between the coefficients in this quadratic form and those appearing in fyy_ (by, ---, by). 

The number of such coefficients is r(r + 1)/2, a number independent of M and N. Con- 
Sequently, if r = 5, or even 10, we have a feasible procedure which is superior to the direct 
method outlined above. Details of this approach and many further results will be found in a 
forthcoming thesis by R. Beckwith. 
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CONSTRAINTS 





The introduction of constraints leads to complications in this case, since the admissible varial 
j 1 
range of the y,, will depend upon the xj. Consider the special case where there are only a small \ns 
number of constraints — for example, one. Let this be feasil 
- 89 
(82) 0 < yin < Xin: (89) 
Since the value of x, itself changes as the y,_4 are determined, this becomes an essential 
state variable. It follows that having transformed the variational problem into the form 
appearing in (78), we now introduce the variable z = Xj,, and define the sequence (90) 
N-1 
(83) fyy(by,b9, ---,by3z) = Max g(bj + ...,b,+ -..) —A 2 bvwd , over | 
os elimi 
monot 
The recurrence relation is then expec 
(84) fy(by, ---,bp;z) = Max [— Ah(yg) + fy_4(by + bygyg,---»2 + Ly(yo) | , 
7 succ 
0 
where L4(¥o) is a linear function in the components of Yo: obtained from the relation %.4* 
Ax, + Yn» and the maximum is over all yg satisfying the constraints independent of xg, and nonlir 
(85) 0< ¥49 <2. (91) 
Let 
SUCCESSIVE APPROXIMATIONS —I 
Let us now consider the case where the criterion function measuring the value of the (92) 
state at time N depends upon more than just a few of the components of Xn To treat this 
case, we wish to use successive approximations. where 
0 
Let } vn be an initial set of policy vectors and \x° ' the corresponding set of state be fou 
vectors. Then 
N-1 POLY 
0 0 0 0 
(86) Jy") =y (x ins ---o%pn) —A Dhlyp). 
k=0 
pointe 
Consider the new problem of maximizing wet 
feasit 
N-1 
ae 0 0 
(87) 559) = 6 (Rago ZaneXgay --+0%yqq) — > DM)» 
k=0 polyn 
: : ; Over | 
subject to the linear recurrence relation 
appra 
” t0 apy 
(88) Se AX, + Yy? 
degre 


and constraints on the y, which are independent of the X,,. 
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As we know, this problem can be treated by means of sequences of functions of two 
variables. Let us then assume that we can consider a problem of this nature solved. Let 
yn \ ‘ x ' be the set of policy and state vectors derived in this way. Since } y? t isa 


feasible policy, clearly 


: 0 
(89) Iy(y') > dy ). 
Now consider the problem of maximizing 


— 0 Fe! 
(90) G(X) > Sony Xgnp X4nr °° Xen) >» (¥,) 


over the y,- This is accomplished as shown above. We repeat this process until we have 
eliminated all of the variables with zero superscripts. Proceeding in this way, we obtain 
monotone approximation. The question of convergence, and to what, is, as might be 


expected, nontrivial. 


SUCCESSIVE APPROXIMATIONS — II 
Let us now mention the point that the preceding results can also be utilized to treat 


nonlinear recurrence relations of the form 


(91) X41 7 UXpr Vn) XQ =C- 
0 0 sale ; ; , ; 
Let y,,X, be an initial approximation. Then we consider the new recurrence relation 
00 0 0 
(92) Xne1 = UX Vp) + Ap(X, — Xp) + BYl¥, — Yn) XQ =C; 
, : , (0 0) 
where A, and B,, are obtained by expanding Q(Xp,Y,) about \Xn»Yn( - Further comments may 


be found in [18]. Much remains to be done concerning the convergence of such procedures. 


POLYNOMIAL APPROXIMATION 

Finally, let us mention an approximation method of an entirely different type. As 
pointed out in the preceding pages, the direct tabulation of a function of k variables for 
k > 4, and at the moment even 3, over the points of any grid of reasonable fineness is not 


feasible. 


How then shall we store the function? One answer is provided by the method of 


Polynomial approximation. Suppose that we approximate to a function of one variable, f(u), over 


over an interval 0 < u < 1 by means of a polynomial of degree 10, usually a quite satisfactory 
approximation. This requires eleven coefficients and yields all values in|0,1]|. Similarly, 

{0 approximate to a function of two variables, f(uy,U9), for 0 < Uy,Ug < 1, by a polynomial of 
degree 10 in uy and ug, requires only 66 coefficients; a function of three variables requires of 


the order of magnitude of 250 coefficients for similar approximation; and so on. 
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Whereas the number of grid-points increases exponentially with the dimension, the 


number of coefficients increases at a much slower rate. How to combine this method with 
the functional equation technique is not a small matter. A preliminary discussion is contained 
in [19], and a great deal of work in that direction is currently being carried on by C. Hast- 
ings; see, for example, | 20]. 
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PROBLEMS 

It has been suggested that this journal might serve as a medium for publishing problems 
in the area of Logistics — with the idea that other persons might be interested in those problems 
and might submit comments. Readers are invited to submit brief statements on applied and 
theoretical problems in Logistics. Address letters to Managing Editor, Naval Research Logis- 
tics Quarterly, Office of Naval Research, Washington 25, D. C. 


24 December 1959 


Dear Sir: 

The letters of McShane and Solomon (NRLQ December 1958) and of Davis (NRLQ June 
1959) commenting on several articles dealing with the Naval Electronics Allocation Problem, 
including Kruskal’s and my article [1]: “The Coefficients in an Allocation Problem” (NRLQ 
June 1958), have recently come to my attention. Most of the points raised in these letters were 
covered in Kruskal’s and my article [2]: “Assigning Quantitative Values to Qualitative Factors 
in the Naval Electronics Problem” (NRLQ March 1959). Still, because of the importance of the 
questions raised, I think it is worthwhile to go over some of these points in the less technical 
context of a letter (especially as Davis’s letter appeared after [2]). 

(1) [1] is not an essay on the Electronics Allocation Problem, but a description of a 
general approach to subjective allocation problems, i.e., to allocation problems in which the 
values of the individual assignments are not given numerically at the outset, but must be deduced 
from “qualitative” information (such as preferences; cf. (9) below). We referred to a greatly 
simplified version of the Electronics Allocation Problem in [1], to illustrate the general approach 
without getting lost in a sea of complications. The application to the real-life Electronics Allo- 
cation Problem was described in [2], not in [1]. This accounts for the fact that the terminology 
in [1] was, according to McShane and Solomon, “less clearly defined than it might be,” and that 
“there still remains the need for a sharper definition of the actual problem.” I hope that this 
sharper definition has been provided in [2]. 

(2) McShane and Solomon state that the MIP “is the basic source of information in the 
problem.” It is true that the MIP is a basic source of information in the problem, and that it is 
the only basic source readily available in organized form, but the whole complexity of the prob- 
lem lies in the fact that there are two other basic factors that must be taken into account, and 
that interact with each other and with the MIP. These factors (goodness of model available for 


installation and of model already installed) are not as readily available in organized form as the 
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MIP, but that does not mean that we can ignore them (cf. pp. 2-3 of [2]). This fact was recog- 
nized when, during the experimental phase of this work, an officer was appointed to decide “tokey’ 
questions in which the various factors indicated conflicting answers. If the MIP were the only 
basic source of information, these questions could have been decided trivially. 

(3) Neither of the two interpretations of the MIP suggested by McShane and Solomon is 
the one used by us. Interpretation (a), that “all priority 1 requirements be fulfilled before pri- 
ority 2”, no matter what equipments are installed and no matter how suitable the available mode] 
is for the two positions in question, would be too drastic an oversimplification. Interpretation (b), 
that total “effectiveness” should be maximized, is too vague to be useful. Effectiveness is a 
derived concept, constructed with the aid of the MIP and other information, and cannot be used 
to give the intent of the MIP. The interpretation used by us is precisely given on p. 5 of [2]. 


It says that all priority 1 positions should be filled before any priority 2 positions provided all 








other factors are equal. 








It would be difficult to deny that the MIP says at least this. That it does not say more 
can be seen from the fact that duly constituted Naval authority has in the past decided “token” 
problems against the MIP, on the basis of the other factors. Such a procedure can also be easily 
justified on a common-sense basis: the MIP contains hundreds (if not thousands) of items, all in 
rank order; it seems intuitively clear that there are some pairs of items that are close to each 
other on the MIP, that “differ very little” in importance. Now suppose that by and bs are posi- 
tions corresponding to MIP items of this kind, by being preferred to by. Let us say that there 
is a fairly good equipment installed in b, that is doing the required job, maybe not in the best 
possible way, but fairly well; whereas in bo there is installed a very poor equipment or no equip- 
ment at all. Suppose moreover that the only equipment currently available for installation is 
eminently suitable for installation in bg, but comparatively unsuitable for installation in by, say 
because of weight, size, or similar considerations. I think it is fairly clear that the available 
equipment should be installed in bg in spite of the MIP. Now the questions arise, how “near” do 
by and bg have to be in order for this reasoning to be valid? Conversely, how “great” does the 
difference between the already installed equipments (and/or between the suitabilities of the 
available equipment to the two positions) have to be to justify going against the MIP? This is 
the crux of the Naval Electronics Allocation Problem. Simply to say that we will follow the MIP 
regardless of the other factors is to beg the principal question. 

Of course in a given practical allocation problem, all other factors are rarely equal. 
Nevertheless the MIP, when interpreted as above, is one of several powerful tools that together 
can be used to arrive at practical allocation plans. I don’t think the MIP was ever envisaged as 


a document that would immediately and single-handedly solve allocation problems. 


(4) Questions asked of the “Board” (the individual or group authorized to arrive at allo- 
cation decisions) involve realistic situations only ( [1], p. 122). Nobody would consider using a 
VHF transmitter to fulfill the requirement of a surface search radar; questions involving such 


a possibility would not occur. 
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(5) We used the word “utility,” in [1], in the following sense: A utility on a set X of out- 
comes is a real valued function v on X, with the property that for all outcomes x and y, x is pre- 
ferred to y implies that v(x)>v(y), and x and y are indifferent implies that v(x) = v(y) (probability 
combinations are irrelevant in this context). Under this definition, I think it justified in a busi- 
ness problem to call dollar profit a “natural” or “physical” utility. If this definition is rejected, 
the word “objective function” may be substituted for “utility.” 

(6) McShane and Solomon make the point that some allocation problems do not possess 
the complicating features with which we deal. This is all to the good. If an equipment is totally 
unsuitable for all requirements and ship classes but one and is the only equipment available, and 
if all positions of this requirement and ship class have voids currently installed, then or course 
there is no problem; one simply installs the equipments on the ships until one runs out of them. 
There are, however, cases in which things are not so simple; and it is to these cases that our 
work has been directed. 

(7) Intransitivity would yield inconsistent numerical results, and thus would be sufficient 
cause to consider our technique inapplicable ([1], pp. 120-121). We did not run across any 
intransitivities in the experimental phase of our work. I think the intuitive meaning of “prefer- 
ence” precludes intransitivities. 

(8) The problem raised by Davis, that the second radar installed on a ship should not have 
the same value as the first, was discussed under the heading “Multiple Allowances” on p. 12 of 
[2]. The conclusion we .eached there may be expressed in Davis’s notation by saying that 
(2)a;; means “the installation of two radar sets on two distinct submarines of the same priority 
class,” and not “the installation of two radar sets on a single ship.” Thus (in Davis’s words) 
“we might hope to rule out diminishing marginal utilities.” Davis’s conclu ‘on that his plan e 
could never be preferred to his plan d is not justified, because he has faile .» take into account 
factors other than the MIP. 

(9) I conclude with a rough explanation of the “logical underpinnings” of our technique. 
Like Davis’s, our approach uses only orderings. Let us concentrate our attention on a fixed 
String of requirements (a “string” is a set of related requirements that for the purposes of 
allocation are considered together; cf. [2], p. 3). Basically, we know only that there exists a 
partial preference order on the set X of all allocation plans for the given string; all this says 
is that it is sometimes meaningful to say that one plan is preferred to another. Our task is to 
find a maximal feasible allocation plan, i.e., a feasible plan such that no other feasible plan is 
preferred to it. (If there is a maximum feasible plan, i.e., a feasible plan preferred or indiffer- 
ent to all other feasible plans, then every maximal plan will also be a maximum plan; but unless 
the given preference order is total, a maximum plan may not exist.) 

We now make certain assumptions about the preference order; the significant ones are 
transitivity and an additivity assumption, which says approximately that if plan A is preferred 
to plan B, and if we add the same assignment a to both A and B, then Aja is still preferred to 


B+a. (This additivity assumption is closely related to the idea of constant marginal utility.) It 
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can be proved that there then exists a utility v on X (cf. (5) above), with the “linear” property 
that v(A+B) = v(A)+v(B). Maximization of v over the set of all feasible plans yields a feasible 
plan that is maximal in the given preference order. Thus if we can succeed in ascertaining v, 
we have reduced our problem to a standard linear programming problem. 

The trouble is that the preference order exists only in the mind of the “Board,” and even 
there it exists only in disconnected pieces. To ascertain the preference order explicitly is out 
of the question. Even ascertaining the utility v precisely would usually involve asking the board 
thousands of questions about the preference order, and would be impossible from the practical 
viewpoint. On the other hand, it is possible to determine v approximately; this is accomplished 
by using the information contained in the MIP and in the answers to a limited number of “token” 
questions. If w is an approximation to v, then maximizing w will not necessarily yield a maxi- 
mal plan, but it will yield a plan that is in a certain sense “close to maximal;” without going 
into details, we may say roughly that the better the approximation of w to v, the closer the 
resulting plan will be to maximality. In any case, the plan obtained by maximizing w will be 
“best possible” with the given information. 

It is important to realize that neither the statement of the problem nor its final solution 
involve the concepts of utility or “military worth.” These are derived concepts, used as com- 
putational tools in arriving at allocation plans that are maximal or near maximal. Certainly they 
can have no direct real life interpretation such as “range” or “search radius.” 

I would like to take this opportunity to thank McShane, Solomon, and Davis for their 
interesting and stimulating comments on this important problem. 

Sincerely yours, 

/sgd/ ROBERT J. AUMANN 


The Hebrew University 
Jerusalem, Israel 
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A COMPREHENSIVE BIBLIOGRAPHY ON OPERATIONS RESEARCH: 1957, Operations 
Research Group, Case Institute of Technology. 123 pp. $2.00. 

Preparation of this bibliography was made possible by a grant from the Office of Scien- 
tific Information of the National Science Foundation. Under this grant bibliographies for 1958 
and 1959 also are being prepared. 

Admittedly this is not a completely perfect bibliography, but despite its imperfections 
it will be of considerable use to anyone working in the field of operations research. 

This bibliography may be obtained from: Mrs. Grace White, Secretary, Case Institute 
of Technology, University Circle, Cleveland 6, Ohio. 

(Reviewed by J. Campbell) 


PORTFOLIO SELECTION, by Harry M. Markowitz. Monograph 16, Cowles Foundation 
for Research in Economics at Yale University. John Wiley and Sons, Inc., New York, 1959. 
viii + 344 pp. 

Professional economists have long since come to expect from the monographs and other 
publications of the Cowles Commission that they will represent the applications of careful 
reasoning to significant problems. The present monograph, which represents the summation 
of Dr. Markowitz's research on the question of efficient diversification of investments, of 
which some has previously appeared in various journals, is no exception to this rule and 
testifies to the careful thought which has so evidently gone into its preparation. 

The monograph is divided into four parts. In Part I Dr. Markowitz introduces the 
general problem of portfolio selections and indicates the nature of his proposed solution. 
Portfolios are rated according to two criteria, expected returue and expected variability of 
these returns as measured by their standard deviation or some related concept. The evaluation 
method is as follows. Ratings on individual securities, both as to expected returns and the 
expected variation of the returns, are obtained from experienced investment councilors. The 
portfolios obtained by combining these securities are compared, and those portfolios which 
are "inefficient''—that is, have both lower expected returns and higher expected variance than 
Some other portfolio in the group—are discarded. From the remaining set of "efficient" 
portfolios the investment advisor then selects that particular portfolio whose combination of 
expected returns and expected variability of return is most suited to the individual preference 


pattern of the investor for whom it is prepared. 
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The remainder of the book is an expansion of this central theme. In the second section 
the author considers the statistical concepts involved in his analysis: averages, expected 
values, variance, co-variance, standard deviation, and the effect of the number of securities 
in the portfolio and of the time period involved. The third section is concerned with the 
question of efficient portfolios and presents a discussion of the geometric concepts involved, 
a method for deriving efficient portfolios, and a discussion of semivariance as a possible 
alternative measure of risk. 

In Section IV, Dr. Markowitz introduces the notion of expected utility and utility 
functions. After a rather lengthy discussion of the pros and cons of the expected utility 
hypothesis, the treatment of utility maximization over time, and the nature of probability 
beliefs, he shows how it may be possible to derive utility functions which adequately represent 
the utility of wealth for investors and to utilize these functions in the selection of portfolios. 
Three appendices give, respectively, the mathematics involved in the computation of efficient 
sets of portfolios, a variant of the simplex method for machine computation of these sets, and 
a discussion of the various axiom systems involved in the expected utility hypothesis. 

Certain minor objections may occur to the reader. In an obvious effort to make the 
volume self-contained the author has devoted considerable space to a rather exhaustive 
discussion of certain elementary statistical concepts, the knowledge of which might well be 
assumed, at least, by persons who would have any serious interest in the work. Ina related 
context there is some doubt in the reviewer's mind as to exactly whom the author considers 
as his reading public, since for the mathematical economist much of the material presented 
would be too elementary or "old hat,'' while, if the book is intended as a practical guide to 
help investment advisors actually select portfolios, some of the mathematical niceties might 
well be dispensed with. One can only conclude that Professor Markowitz is addressing 
himself to both groups at once—a difficult task and one which he has done as well as could be 
expected. The whole relevance of the utility discussion might be questioned if the problem is 
one of selecting a portfolio for an organization, certainly one of the more common problems 
which is encountered. 

One final point. Although Professor Markowitz confines himself solely to the problems 
of security and portfolio analysis, his method has certain possibilities for use in similar 
decision problems — selection of an optimal weapons system or set of development projects, 


for example — and may well be valuable to persons involved in these fields. 


(Reviewed by W. G. Mellon) 


MATHEMATICAL MODELS OF OPERATIONS RESEARCH, by Dr. Thomas L. Saaty. 
McGraw-Hill Book Company. $10.00. 

This book by Dr. Saaty fills a definite gap in the existing Operations Research litera- 
ture and must be considered a major contribution. A surprising amount of valuable informa- 


tion has been gathered together in this single volume. 
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To fully comprehend Dr. Saaty's book a background in calculus beyond the usual first 



































course is necessary, as is a minimum knowledge of matrix theory. This is not to say that the 
es less sophisticated mathematician will not benefit from it, nor does it indicate that the more 
sophisticated mathematician will not find the book both informative and thought provoking. 

\ Examples are used throughout the book, and problems are provided. The book also 
contains an excellent list of references. As a pioneering effort, ''Mathematical Methods of 
Operations Research" will doubtless aid future authors and will influence their work in this 
field. This volume belongs in the library of anyone interested in the science of Operations 
Research or the techniques it employs. 


It is interesting to know that this volume is being translated into the Japanese language. 


ent (Reviewed by G. G. O'Brien) 
ent 
_ METHODS OF CORRELATION AND REGRESSION ANALYSIS, by Mordecai Ezekiel and 


Karl A. Fox. John Wiley and Sons, New York, 1949. xv + 548 pp. 
, The whole tone of this book, the third edition of Ezekiel's well known Methods of 


Correlation Analysis, is set by the first paragraph of the preface: "'Thirty years have elapsed 





since the original edition of this book was written...(but)...the basic elements of correlation 
ed analysis continue unchanged.'' Just how ridiculously untrue this is may be seen from the 


contents of two recent issues in the Wiley Publications in Statistics series, E. J. Williams' 


we 


d Regression Analysisand F. §. Acton's Analysis of Straight Line Data. 








The book is entirely nonmathematical and is obviously aimed at a group of research 
ht workers (should they exist) who want to know what to do but not why they doit. The additions 
to the original edition add little to its usefulness, and the chapter on error formulas in time 
be series is definitely misleading, as it almost entirely ignores the important results of the last 
fifteen years. The book remains valuable for the number of worked examples, but it is 
doubtful if anyone owning the earlier edition will benefit from obtaining this new edition. 


(Reviewed by C. W. J. Granger) 


ems 

ANALYSIS OF STRAIGHT-LINE DATA, by Forman S. Acton. John Wiley and Sons, 
” New York, 1959. xiii + 267 pp. 

REGRESSION ANALYSIS, by E. J. Williams, John Wiley and Sons, New York, 1959. 
on) ix + 214 pp. 

These two volumes belong to the fast-expanding and mostly excellent Wiley Publications 
y: in Statistics and are both well up to the average standard for the series. Although the subjects 

dealt with in the two books often overlap, their individual methods of presentation, together 

ais with the numerous topics not in common, make each worthwhile. 
1a- Professor Acton has produced an important book of by no means trivial contents, even 


though it is hidden behind a trivial-sounding title. The difficult task of introducing and 





differentiating between all the various regression situations is carried out most successfully 
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and carefully. The book is aimed at intelligent research workers with a fair amount of 
mathematical training, and any such worker who is likely to be meeting the problem of fitting 
a line to data often in his life will find reading this book extremely worthwhile. 

The approach throughout is entirely honest, the assumptions being required are care- 
fully explained and the limitations of any method admitted. The philosophy of an approach is 
always given due attention, and fitting a line blindly by least squares is never advocated. 
Nonparametric and low-arithmetic methods are always prominent, and in these fields, together 
with the idea of transformation of variables, several interesting but previously unpublished 
results due to Professor Tukey are introduced. Two unusual and noteworthy aspects are the 
consistently good and useful diagrams to be found throughout the book and the two valuable and 
carefully prepared supplements, one on finite populations and the other concerned with fitting 
orthogonal polynomials for unequally spaced points. 

The first half of Acton's book, which deals with the case of one of the variables known 
without error, is considerably more valuable than the second half, which deals with a miscel- 
lany of other cases. The only criticisms that can be levied, apart from the unevenness of the 
last few short chapters, are the lack of numerical examples with which a reader may try out his 
new skills and the surprising omission of a discussion of the idea of functional relationship, the 
paper by Lindley in Biometrika (1953) not even appearing in the supplementary bibliography. 

Both of the books under review are written against a background of wide practical 
experience, and neither could be strictly called a mathematical work, as few of the statements 
are proven. This lack of mathematics makes Dr. Williams’ task more difficult than that 
tackled by Professor Acton, as he is dealing with intrinsically more complicated problems. 

For the first half of his book, Dr. Williams introduces all the classical regression 
problems with clarity but generally with a less ''modern" approach than that used by Professor 
Acton. However, in the second half, the book proves its true importance by a series of well- 
done chapters on topics not usually found in a textbook. Here the treatment of heterogeneous 
data, simultaneous regression equations, discriminant function, and functional relations are 
all considered . 

It is felt that each of these two books dealing with different aspects of an oft- 

misrepresented subject will prove to be very useful to experimenters, practical statisti- 
cians, and teachers of statistics. 


(Reviewed by S. N. Afriat and C. W. J. Granger) 
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NOTES ON OPERATIONS RESEARCH 1959, assembled by the Operations Research 
Center, M.I.T. The Technology Press, Cambridge 39, Massachusetts. 256 pp. $4.00. 

The eleven chapters of these notes, assembled by the Operations Research Center, M.I.T., 
cover probability, search, Markov processes, queuing, control processes, organization of 
operations research groups, sequential decision processes (the most enjoyable chapter which 
includes some novel features), reliability and maintenance, information theory, production 
scheduling, and simulation of random processes. 

From the preface one gathers that some of the ideas were used as lectures to persons 
from NATO countries assembled in Brussels in August 1959. Most of the material presented 
is found in existing operations research literature, a few examples of which are given in the 
notes. It is asserted that no pretense is made to cover all aspects of operations research, etc. 
The notes may be of interest to those inclined toward a quick familiarity with methodology. A 
chapter which developed problems from one operation and showed solutions by use of these 
methods would have been highly desirable. This would have improved the apparently fragmentary 
structure of these notes (and the included problems) and given the reader a unified outlook as 
to their utility. 

(Reviewed by Thomas L. Saaty) 


LOGISTICAL SUPPORT OF THE ARMIES, Vol. Il, by Roland G. Ruppenthal. U. S. 
Government Printing Office, 1959. 511 pp. plus glossary and index. $4.50. (This is the forty- 
third volume published in the series, UNITED STATES ARMY IN WORLD WAR II, and the 
fifth volume in the subseries, “The European Theater of Operations.”) 

This admirably organized work emphasizes the limitations logistics imposes on strategy 
and tactics through the central problems of port rehabilitation and management, movement con- 
trol and transportation, supply of rations and fuel, equipment, ammunition and manpower. 
Throughout, the treatment is clear, objective, and analytical. 

As in Volume I, the author cites repeated instances of friction between tactical and logis- 
tic commanders. While some of these were undoubtedly due to personality conflicts, most 
evolved from the basic problem of providing a flexible organization to deal with a rapidly devel- 
oping situation. 

Starting with Chapter II and throughout the book there is continuing discussion of the 
difficulty of relating tactical plans to the twin dynamics of port development and transportation. 
Ship berthing and unloading, cargo clearance, depots, motor and rail capacity, tactical move- 
ment, strategic objectives, strategic momentum, availability of shipping, priority of unloading, 
and selective unloading — all these are related and are regenerative. Once a bottleneck has 
been established in port cargo clearance, a vicious spiral of ever-increasing trouble, danger, 
and waste develops. This ever-recurring story began at Noumea in 1942 and was repeated 
throughout the whole war. Its basic causes and effects will continue in future war wherever the 


movement of forces in combat is required. 
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When one studies the complicated relationship between effective combat power, combat 
unit build-up, logistic support and build-up, transoceanic shipping, port capacity, etc., discussed 
in Chapter V and subsequently, he wonders how thoroughly these lessons have sunk in on those 
who are conducting research in linear programming and war gaming as related to logistic- 
strategic policy, force structure, and budget allocations. The cost of such research and war 
gaming is trivial in relation to the stakes involved and to the cost of weapons research. 

Chapter IX, on ammunition, provides specific illustration of the direct limitation of 
tactical operations by logistic shortages. This account of alleged and real ammunition shortages 
shows not only the great importance of good logistic planning factors but also the vital impor- 
tance of fire discipline and of logistic discipline. 

This and the description of the problems arising from the inability to forecast attrition 
factors in equipment and tanks make an interesting contrast to Marshal Rommel’s accounts of 
his African Campaigns. What was considered a scarcity by the Americans would have been an 
extraordinary luxury to the German leader. 

The discussion of manpower in Chapter XI brings out the vital importance of quality and 
versatility of personnel as being the basis for flexibility to meet new or changing conditions. 
And, when one rereads the chapter, it becomes obvious that enormous administrative and log- 
istic effort was misdirected and misused. In other words, the whole complex of the logistic 
administrative snowball and its relation to the eternal problem of centralization versus decen- 
tralization of authority is exposed. Disagreement as to the basis for plans caused trouble. 
Ruppenthal points out: “Unfortunately the manpower problem, like the ammunition problem, was 
needlessly aggravated by the lack of mutually understood rules of accountability.” 

Chapter XIII, “Tactical Logistical Organizational Aspects of the Last Offensive,” gives 
an excellent appreciation of the value of alternate plans and of the manner in which logistic 
capabilities and flexibility enable a strategic exploitation of a tactical success. “Success in 
supporting the final drive can be attributed in large measure to the flexible plans which COMZ 
organizations had worked out for the expansion of transport in the forward area and the rapidity 
with which rail bridges were installed. The decision to prepare simultaneously for alternate 
offensives with maximum strength in both north and south proved wise, for these preparations 
made possible the prompt exploitation at Remagen and adequate support of the crossings in the 
Frankfurt-Mannheim sector.” 

Throughout the book Ruppenthal repeatedly points out how the clear resolution of key 
problems of command and organization was evaded or postponed until after serious damage was 
done. His comment, “A genuine theater headquarters eventually emerged in Europe, but only 
as the result of the adjustments incident to the end of hostilities,” is an interesting contrast to 
the relative clarity and permanence of Admiral Nimitz’s CINCPOA, CINCPAC commands of 
1943-45. 

This problem of command and organization will become particularly important tomorrow 
when vital command decisions will have to be made much faster because of the great increase in 


speed and power of modern weapons. This great increase in speed of events and decision will 
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demand such rapid collection evaluation and display of essential information and intelligence that 
Ssed special equipment must be designed and installed long before events require its use in crisis or 
combat. 

Since normal lead time in the development of such equipment is several years, it is highly 
doubtful if major changes in command responsibility can effectively be made after a crisis devel- 
ops. It is, therefore, vital that the logistic aspects of future combat operations be thoroughly 
analyzed now in order that the command responsibilities can be carried out. The basic principle 
Ages can be stated quite simply: In modern conflict the man who understands and controls the capa- 

- pilities and location of electronic command equipment can in effect decide who will wield actual 


command authority. 


n Chapter XVI, “Supply in the Last Months,” and Chapter XVIII, “The End of the Replacement 
of Problem,” vividly bring out the principle that initial underplanning results in overplanning and 
in the eventual accumulation of great and unnecessary surpluses. For example, as late as 25 April 


ETOUSA cabled that it wanted the flow of men from the zone of interior maintained. Ruppenthal 
and says, “At the very moment that ETOUSA was painting this very gloomy picture, it was accumu- 
lating replacements far in excess of its immediate needs, and actually faced a problem of how 

- to use the surpluses.” 

The lesson is clear; once mutual confidence between superior and subordinate is destroyed 
n- the pendulum can make wide swings between acute shortage, true privation, and reckless over- 


estimates and wastage. 





We re eS Se St Sh a A ae Se 
This reviewer feels that in analyzing the logistic problems presented by excessive or 
'S ill-timed combat force build-up and by other decisions aimed at pressing home offensive 
operations as hard and as rapidly as possible, one must be careful to keep the true logistic 
objective in mind. The true logistic objective always should be “to attain maximum sustained 
1Z combat effectiveness.” 
dity With this in mind one can view the exasperations and frustration of logistic commanders 
2 and the disruption of logistic plans with equanimity, provided the actions resulting in such 
ns frustration and disruption actually did contribute to the combat objective. When, however, these 
the disrupting actions detracted from the attainment of the combat objective we then have something 
else entirely. 
It requires very careful analysis and mature judgment to discern the difference when one 
was is making command decisions in these areas. This is precisely the reason “operators” must 
ly study logistics and logisticians must study strategy and tactics and have the “Perspective of 
to Combat Command.” 
The concluding chapter, “In Retrospect,” which should be required reading for all officers 
‘row aspiring to high command, brings to a fitting and perceptive close the best organized and most 
se in 


analytical discussion of World War II logistics yet published. In this conclusion Ruppenthal 


clearly sums up the most important facets of the logistical problem in modern war, specifically: 
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Its vital influence upon strategy and tactics; 


The troubles which inevitably ensue from uncertain or faulty 
command structure; 


The size and complexity of modern logistics; 
The principles of logistic limitations; 
The problem of balance of forces; 


The principle of logistic reverberation — the “Dynamics of 
Logistics;” and, finally and very appropriately 


The importance of assigning officers of proven character and 
ability to logistics tasks. 


As we look into the future with all its technological advance, there seems nothing to indi- 
cate that the importance of these fundamentals, nor the need for military commanders to think 
carefully about them, will in any way diminish. On the contrary, the congressional hearings and 
public debate on defense issues in the winter of 1960, with their emphasis on the enormous and 
complex logistic effort involved in “Airborne Alerts,” “Hardened Missile Bases,” and “Defense 


Highway Bridge Heights,” indicate that logistics considerations will dominate the future. 


(Reviewed by Henry E. Eccles) 


PRINCIPLES AND APPLICATIONS OF RANDOM NOISE THEORY, by Julius S. Bendat. 
John Wiley and Sons, New York, 1958. xv + 431. 

INFORMATION THEORY AND STATISTICS, by Solomon Kullback. John Wiley and Sons, 
New York, 1959. xvii + 395. 


Statistics has often benefited from considerations of problems occurring in other sciences, 


but rarely has the combination been so fruitful to both subjects as the partnership between con- 
munication engineering and statistics. This partnership has produced not one important field, 
but two — the spectral approach to time series and information theory. Both of these fields may 
be thought of as having been developed since the war and they have matured sufficiently to merit 
expository textbooks of a fair size. 

The two books under review are primarily directed at different audiences, Bendat’s book 
being essentially for communication engineers and Kullback’s for statisticians, but members of 
each field will benefit from both. 

Bendat’s book is sound, informative and purposeful, being careful but not over-fussy with 
mathematical rigour (a type of book, incidentally, which is often found in the applied sciences 
but which is rarely available to the less fortunate mathematician or statistician). 

The concept of a spectrum is carefully brought out in the first few chapters; the combined 
notions of ensemble averages (over several samples from the same time series), averages ove! 
time, ergodicity and stationarity being better combined here probably than in any other book. 0! 
the other hand, the difficulties of estimating the spectra, which has produced a flurry of public* 
tions recently (see (1), (2)) is entirely ignored, although a whole chapter is devoted to statistical 


errors in autocorrelation measurements. 
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Having introduced the most important tools of the field, i.e., the power spectrum and the 
filter, the author then considers a variety of more or less unconnected topics such as analog 
computer techniques of measuring the root mean square ensemble output errors of (possibly 
nonstationary) random inputs characterized by exponential-cosine autocorrelation functions, 
envelope detection and the zero-crossing problem. These sections contain several new devel- 
opments and derivations together with suggested applications. 

Statisticians will find a most refreshing and important feature of the book the fact that 
nonstationary series are discussed almost as frequently as the more common stationary series. 
In particular, the chapterson “Optimum Linear Prediction and Filtering” and “Optimum Time- 
Variable Filters” present the techniques used by communications engineers on nonstationary 
data in an easily assimilated form and various new results by the author are included. The 
contents of this book bring out clearly how the research concerning time series has split again 
into two parts: communication engineers being primarily concerned with filtering and non- 
stationarity, while statisticians have been turning their attention mostly to methods of estimating 
and interpreting spectra and dealing with multiple time series. 

Such a split has also occurred to some extent in the field of information theory, although 
here the statistical results are moving ahead faster. 

Kullback has written a book aimed directly at the professional mathematical statistician, 
introducing and examining the properties of a measure of the amount of information contained 
ina set of data and then discussing various of its cases in the field of statistics. 

In the first five chapters of the book, the author gives to information theory the same 
background of sound mathematical rigour that Doob, in another volume in the same series, has 
given to stochastic processes. The algebraic properties, inequalities, limiting properties and 
estimates of the information function as introduced by Shannon are fully investigated but there 
is no discussion of the various generalizations of measurement of information that are now 
available (e.g., see (3), (4)). 

The applications start with chapters on tests of hypotheses involving multinomial and 
Poisson populations and these are followed by an excellent and important chapter on contingency 
tables. The application of the idea of information to such a well-tried method and the interpre- 
tations and extensions that result clearly show the potency of the theory. The later chapters 
deal with tests concerning the multivariate normal populations and multivariate analysis with 
various hypotheses. 

Although the fields with which the two books are concerned may be said to have a common 
ancestry, it is perhaps surprising how studiously they ignore one another. Bendat never men- 
tions the word “information” and Kullback mentions “stochastic process” once in his introduc- 
tion but never again. Since time series and the analysis of noise have reached a stage in which 
tests of hypotheses are under consideration, e.g., stationarity, smoothness of a spectrum, inter- 
dependence of two series, the present would seem to be the opportune time to apply the results 
of information theory. 

Little of use to the time series analyst has yet been produced, but some of the possibilities 


of the field are brought out in a recent paper by the Russians I. M. Gel’fand and A. M. Yaglom (4). 
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It has long been recognized that the degree of dependence between two series is measured at a 


given frequency A, by 
2 2 
22> , 0< CAd< 1 
f(A) g(a) 


where f(a), g(A) are the spectra for the two series, c(A) is the co-spectrum and q(d) the quadra- 
ture spectrum. Gel’fand and Yaglom have put this knowledge on a better footing by showing that 
the average amount of information per unit of time concerning one series derivable from the 


other is given by 


oo 
i(x,y) = -5 | 6 [1 . C(a)| dv 

or, at frequency 
ix, ¥;A) = 5 log [1 - C(a)] ; 

For samples of length n from each of the two time series the natural estimate of 

i(x,y;A) is 
i(x,y;A) = -5 log [1 - C(a)] 

where 


a aa a2 
&) = © (a) +q“(a) 
f(A) g(A)” 
the cap here indicating an estimate. If the estimates used are similar to the Tukey estimates, 


N. R. Goodman (5) has shown that under certain unrestrictive conditions €(a) has frequency 


function 





n = co 2 
ans > ‘(nk eek) (1-22 
r(n)r(n-1) a0 r“(k+1) 


where r = 1 - C()). 
Thus if the two series were independent at (or near) the frequency \, C(a) has frequency 


function 
p(z) = 2(n-1) 2(1-z2)""2 : os < i 
and so confidence limits may easily be given to i(x,y;A). 
Similarly, defining 


m 
i(x,y) = 5 rm log [1-C(a4)| » Ai41 - Ay = Const, all i 
i=l 


and if the smoothing devices used on the spectral estimates were perfect and the series truly 
independent, i(x,y) would be the sum of m identically distributed, independent variables and so 


approximate confidence limits may be obtained. Thus, using these results, a logical measure 
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’ estimate established. Should the null hypothesis of independence be rejected, the most likely 


frequencies can be tested. 
The ease with which such a result has been found would suggest that further applications 
of information theory to time series analysis ought to produce important results, particularly 


in the testing of hypotheses. 
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Reviewed by C. W. J. Granger) 


NONLINEAR PROBLEMS IN RANDOM THEORY, by Norbert Wiener. John Wiley and 
Sons, New York, 1958. 

This book contains the plan of an instructive and worthwhile monograph about a tool 
invented a long time ago by Professor Wiener. The book is intended to enlighten and interest 
electrical engineers and certain natural scientists whose investigations include random stimu- 
lation of nonlinear systems. The plan has been crudely executed; apparently Wiener undertook 
hd serious proof reading or reworking of the notes. Important gaps, particularly in the basic 
lirst four lectures, will make the going difficult for those with no previous familiarity with the 
material. Nevertheless, at half the price this monograph would be a welcome indication of con- 
tributions long available only in mathematical journals and laboratory memoranda. 

The first four lectures (Chapters) develop the methods. Lecture 1 introduces The Brown- 
ian motion (bm) random process on the unit interval and the corresponding stochastic integral 
and some of its properties. Lecture 2 develops some properties of “homogeneous polynomial 
lunctionals (of the bm process) and their averages.” The lecture concludes by noting the exten- 
sions of the bm process to the infinite interval and complex values. 

The next two lectures are devoted to the representation of arbitrary integrable functions 
by expansions in increasing orders of orthogonal functionals of the bm process. The autocorre- 


lation function is given for this representation and its applicability to questions of spectra for 
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random linear and homogeneous quadratic phase modulation noted. This line of thought is cop- 


bad 


PONE; 


tinued through lectures 5, 6, and 7, in which some computations are undertaken for the two cage; to th 
mentioned. Lecture 8 contains some interesting remarks about sets of coupled nonlinear oscij- 4 Wier 
lators. It also contains some analysis of a particular kind of quadratic phase modulation by 1942 
noise, leading to a spectrum said to resemble qualitatively some spectra actually measured. now | 
No attempt is made to relate the model (equations 8.5 and 8.6) to any physical mechanism. cult | 
Lectures 10 and 11 are devoted to the proposition that bm inputs and representation by expan- when 
sions in ascending orders of polynomial functionals are appropriate for the analysis and synthe- filter 

of Pr 


sis of nonlinear systems. The other five lectures contain miscellaneous comments On quantum 
theory, coding, and statistical mechanics. 

As an expositor, Wiener sometimes presents an infuriating contrast between loquacious- 
ness anent trivialities and silence on important points. In the first lecture (page 10) over a pag 
is devoted to some elementary calculus manipulations, while the immediately preceding switch 
from an average over the space of random functions to an average over the range of the random 
process is accorded three lines. Though the correspondence between measure in the function 
space and the probability density for a given Brownian path is implicit in the construction, no- 
where is it actually mentioned. 

On pages 21-22 of Lecture 2 is introduced an important correspondence, that between a 
transformation of the functions appearing as the kernals of the bm functionals and a measure 
preserving transformation in the space of bm random functions. The ideas at this point appear 
to have gone through a random, nonlinear mutilation. A non-initiate is unlikely to make much 
of this. A similar muddle occurs on pages 40-41, where “closure” of the set of functionals ofa 
given order is in question. 

At certain points things which seem clear to the lecturer are not likely to seem so to the 
average target of the monograph. A trivial instance occurs at the start of Lecture 2, where 
some of the properties asserted to follow from Lecture 1 do only in the sense of being true of 
the objects defined in that lecture rather than because of any explanatory material. In fact, an 
opportunity is lost here to point out that the properties in question are properties of a set of 
normal variates with zero means. A less trivial case arises in Lecture 4, page 43, in connec- 
tion with the previously mentioned relationship between transformations of the functions appeal 
ing as the kernals of functionals and measure preserving transformations in the bm function 


space. Wiener’s “I think that it is clear .. .” or “I think that it is also clear . . .” are unlikely 





to comfort the unenlightened reader who may feel that anything clear should be capable of a 
clear explanation and who would prefer one to being fobbed off by such remarks. In such a 
monograph as this, one expects informal but informative comment. Some is provided — the 
remarks about suitability of the tools for nonlinear system analysis (pp. 27, 37, 88-90) are 
illuminating — but not as much as the informal style might have led one to hope for. 

The equations contain many obvious errors, most of which, one feels, were present in 
the manuscript and are not due to the printers. 
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A final comment about actual application of the tools. The only application carried out is 
to the phase modulation noise spectra considered in Lectures 5 through 8. Having in mind that 
Wiener’s work on the bm process goes back to 1922 or so and that his Rad. Lab. monograph of 
1942 already introduced expansion in functionals of the random process, one had hoped that by 
now he or his students at M.I.T. might have worked by these methods at least one really diffi- 
cult problem of nonlinear system analysis, for example that of finding the amplitude distribution 
when a bm input is filtered, then applied to a linear half-wave rectifier followed by a nonisolated 
filter. Serious progress on some such problem would be the surest demonstration of the value 


of Professor Wiener’s tools, for the tools do seem worth more cultivation. 


(Reviewed by Worthie Doyle) 
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